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Di-neutron clustering and deuteron-like tensor 
correlation in nuclear structure focusing on ^^Li 
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"^ . Abstract "LiisaBorromeannucleus, where two out of three objects as ^Li + « and 

l/^ ' two neutrons independently do not form bound systems. Two neutrons should form 

a di-neutron cluster in the nuclear field generated by the ^Li core nucleus. We treat 
di-neutron clustering by solving the two neutron relative wave function precisely by 
using the bare nucleon-nucleon interaction so that the spatial clustering structure is 
obtained quantitatively within the whole ''Li nucleus. This di-neutron clustering is 
an essential dynamics to form the halo structure by making a compact di-neutron 
^ ' cluster, which distributes loosely around the ^Li core. This concept of di-neutron 

clustering should be clearly distinguished from the BCS pairing correlation, where 

^ no consideration of spatial clustering is made. The di-neutron clustering is a new 

►-j>. , concept and is a general phenomenon in neutron skin and neutron halo nuclei. 

T^ ' This quantitative description of di-neutron clustering has made it necessary to 

f^ . consider another important deuteron-like tensor correlation, which is caused by 

^r ' strong tensor interaction in the nucleon-nucleon interaction. The tensor interaction 

originates from pion exchange and known to provide large attraction to form the 
C^^ ' ^He nucleus. The unique feature of the tensor correlation is to make highly cor- 

^? , related deuteron-like excitation, which interferes with shell model like structure in 

a unique way. This dynamical effect removes the magic number effect and makes 
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easy the participation of the s-wave neutrons. Hence, there are pairing and deuteron- 
hke tensor correlations in addition to the mean field structure in ^Li. The combined 
system of two additional neutrons with the correlated ^Li provides the halo phe- 
nomenon, in which the di-neutron clustering develops with the help of large i-wave 
component caused by the deuteron-like tensor correlation. 

These two effects, the di-neutron clustering and the deuteron-like tensor correla- 
tion, are quite new and essential to provide the halo structure of "Li. In this lecture 
note, we would like to introduce these two new concepts in a systematic manner 
and fill a gap between the halo phenomenon and the microscopic reason for this 
interesting phenomenon. 



1 Unstable nuclei and the halo structure of ^^Li 

We are in the era of being able to study experimentally unstable nuclei up to the drip 
lines and even up to super-heavy nuclei. We are able now to provide precious in- 
formations for astrophysics and cosmo-physics directly from experiment. All these 
activities started with the discovery of the halo structure of "Li H]. Hence, there 
are many experimental data on "Li and surrounding nuclei. First of all, we would 
like to discuss the characteristics of "Li found by various experiments. We should 
present all the existing experimental facts about the Li-isotopes, in particular ^Li 
and i^Li. 

The formation of the halo structure in "Li is a difficult subject to understand 
theoretically. The fact that the binding energy is extremely small urges us to develop 
a theory to handle continuum states as precisely as bound states. This fact of small 
binding energy of the last two neutrons forces us to consider the pairing interaction 
at low density, which leads to the concept of the di-neutron clustering phenomenon. 
Most fascinating physics necessary for the quantitative understanding of the halo 
structure turned out to be the discovery of the deuteron-like correlation caused by 
the strong tensor interaction. Theoretically we have experienced these conceptual 
developments on the theoretical framework to treat the Li isotopes. We discuss first 
the experimental facts on the Li isotopes and then discuss the theoretical tools for 
the understanding of "Li. All these theoretical tools developed for the Li isotopes 
are to be used for many new phenomena found in unstable nuclei. 



1.1 Experimental facts on ^'^'^Liand ^^Li 

There are many experimental data indicating that the "Li nucleus has a halo struc- 
ture. A schematic picture of the halo structure is shown in Fig[T] The halo structure 
of "Li was discovered as an anomalously large reaction cross section of this nu- 
cleus with target nuclei. The matter radii extracted from the reaction cross sections 
and other standard methods are shown in Fig|2] We can see from this figure that 
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Fig. 1 A schematic image of 
halo structure in "Li. 




the matter radius of ''Li is much larger than that of the neighboring nucleus ^Li. 
The matter radius of "Li corresponds to those of medium mass nuclei. If we were 
to pick up nuclei, whose radii are suddenly increased from the neighboring nuclei, 
they are ^He, "Be and '^Be in addition to "Li. Detailed studies on these sudden 
jumps of the matter radii made these phenomena as caused by the halo formation. 



Fig. 2 Observed rms matter 
radii of He, Li and Be iso- 
topes. Shown by the dashed 
curve is the one proportional 
to A^'^ . Those points marked 
by arrows are the ones with 
extraordinary large radii. Data 
are taken from Ref. {!]. 
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There are many experiments performed on the halo structure of unstable nuclei. 
Out of all these experiments, those of "Li are very interesting due to important 
nuclear many body physics behind the experimental data. In order to find the rea- 
son of the sudden increase, experiments were performed for the neutron separation 
momentum distributions. The momentum distributions of neutron separation were 
measured experimentally by Kobayashi et al. i2J by bombarding unstable nuclei on 
some stable target nuclei. The momentum distribution for the case of "Li is much 
narrower than other cases. The narrowness of the momentum distribution is related 
with large extension of the neutron distributions due to the uncertainty principle. 
Hence, this is a data showing directly the halo structure or at least a large spatial 
extension of the neutron distribution in "Li. As for the structure of "Li, there was 
an experiment performed by Simon et al. 13J, who bombarded "Li on a carbon tar- 
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get and measure the momentum distribution of the '"Li fragments. From the shape 
of the distribution, the (lsi/2)^ contribution to the mixture of (1^1/2)^ and (Opin)^ 
components was determined to be (45 ± 10)%. 

If there is a halo structure, we may expect an interesting excitation mode of "Li. 
If the core nucleus ^Li is surrounded by neutrons in the nuclear halo region, there 
is a possibility of making oscillation of the core nucleus in the neutron sea. This 
is called a soft dipole resonance to be excited by photo-disintegration as shown in 
Fig. [3] in |5|. The excitation energy is expected around a few MeV as compared 
to the standard giant El resonance of a few tens MeV. There are several experi- 
ments of Coulomb excitation of "Li. We show the most recent experimental data 
on Coulomb excitation taken by Nakamura et al. in FigUiSl. There is a bump struc- 
ture just above the threshold energy. However, it is still debated if the bump structure 
is caused by the soft dipole resonance or not due to the complicated nuclear structure 
of "Li. There are many other experimental data on "Li as the magnetic moment 
and quadrupole moment. These experimental data will be presented together with 
theoretical results later. 



Fig. 3 A schematic spectram 
of El excitation of soft dipole 
resonance and giant dipole 
resonance modes of halo 
nuclei. Shown above in the 
left is a schematic view of the 
soft dipole resonance (SDR) 
and the one in the right is that 
of the giant dipole resonance 
(GDR). 
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We should also see experimental data for '"Li, which is the neighboring nucleus, 
although this nucleus is unbound. We list the neutron single and double separation 
energies in "^ 10,11 j^j jjj x^^jg i xhe two neutron separation energy in "Li is very 
small as 0.32MeV. The separation energies in ^Li are large and this nucleus should 
be considered as a standard shell model type nucleus. In '"Li, the single neutron 
separation energy is —0.3 MeV, indicating this nucleus is not stable. This resonance 
structure seems to have positive parity and is assigned to have the {p3/2)k{pi/2)v 
structure. At the same time, there are few informations on virtual states. There is 
an experimental data on neutron scattering with ^Li at the threshold energy |7 1. The 
scattering length comes out to be a ~' —20 fm, which is comparable with the one 
of neutron-neutron scattering a = —18.5 ±0.4 fm fS). Hence, the large scattering 
length indicates the existence of strong attraction in the s-wave channel, which is 
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Fig. 4 Coulomb excitation 
of "Li in "Li+Pb at 70 
MeV/nucleon as a function of 
the 'Li+n+n relative energy. 
Details are explained in the 
original papar (6). 
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closed to the condition of forming a bound state in the free space. This means that 
the s-state structure appears close to the threshold energy of ^Li + n. 



Table 1 Single and double neutron separation energies in unit of MeV in ''""Li. 

52„[MeV] 5„[MeV] 



Nucleus 
"Li 
"'Li 
"Li 



0.32 
6.10 



0.62 
-0.3 
4.07 



1.2 Theoretical studies on the halo structure in ^^Li 



The halo structure was completely new in nuclear physics community. Hence, there 
were many theoretical studies to describe this interesting phenomenon. We have 
recognized immediately that the standard shell model approach badly fails due to the 
fact that the two additional neutrons in ''Li ought to enter in the pin neutron orbit 
but not in the sin orbit due to the N=8 magic structure. Hence, most of theoretical 
studies introduce some phenomenology to bring down the si/2 orbit. For example, in 
the work of Thompson and Zhukov ||9|, they treat ^Li as a core and add two neutrons 
by taking state dependent neutron-core interactions. The additional attraction for s- 
wave component makes the (^i^)^ state energetically close to the {pi/2)^ state. In 
this case, the (^1/2)^ state has a large component in the ground state, which provides 
the halo structure for "Li. 
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There is a theoretical study on the pairing property and the E 1 excitation in "Li 
by Esbensen and Bertsch ifTOl . In their study, it is essential to bring down the iw2 
orbit to reproduce the experimental El excitation spectrum. As for the pairing cor- 
relation, there are many studies to describe "Li as the BCS state. In the study of 
Meng and Ring [UTI . they describe "Li in terms of a relativistic Hartree-Bogoliubov 
model. In this study, they can include the continuum effect in their pairing correla- 
tions. In the relativistic Hartree-Bogoliubov model, the s-wave contribution comes 
out to be about a quarter of the p-wave contribution for the paired two neutrons. We 
need more participation of the s-wave component as compared to the finding of the 
experimental data of Simon et al. |[3|. 

There is another interpretation on the halo structure as due to deformation. In the 
work of Varga, Suzuki and Lovas [12|, they try to break the ^Li core and introduce 
the cluster structure. The wave function of "Li is written as ^He+f +4« and take 
the interaction among them by a phenomenological central interaction. In this way, 
they can introduce the effect of the deformation and pairing correlations among the 
nucleons. The deformation effect provides a large matter radius and some s-wave 
component in the wave function. 

The theoretical challenge on the halo structure is therefore summarized as fol- 
lows. There are many indications that the s-wave component is very large in the 
ground state wave function. Hence, we have to find a mechanism to bring down the 
Si/2 orbit with the amount to wash out the N=8 magic structure. The pairing prop- 
erties are also very important to cause admixture of (pi/2)^ and (^i 72)^ states. In 
the halo nucleus, we ought to consider the di-neutron pairing correlation in a small 
nuclear matter density. All these new phenomena should be understood in terms of 
the many body framework with the nucleon-nucleon interaction. 



1.3 Nucleon-nucleon interaction and the deuteron and the 
di-neutron system 

We should learn the properties of the nucleon-nucleon interaction in order to un- 
derstand the halo structure in "Li. To this end, we would like to show the central 
and the tensor interactions in the ^^i channel of the AV8' potential Iil3i . which are 
shown in Fig. |5] In the central interaction, there are strong hard core (short range 
repulsive interaction) and intermediate range attraction of moderate strength. As for 
the tensor interaction, the long range part drops with the pion range, while the short 
range part increases until 0.2 fm and goes to zero at the origin due to the form-factor 
coming from the nucleon finite size. On the other hand, we have the similar structure 
for the central interaction in the '^o channel, where there is a strong hard core due 
to the short range quark dynamics. In this channel, there is no tensor contribution 
due to zero total spin. The deuteron-like tensor correlation is produced by the A^A^ 
interaction in the ^^i channel, while the di-neutron clustering is produced by the 
A^A^ interaction in the '5o channel. 
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Fig. 5 Central and tensor interactions of the AV8' potential [13J in singlet even (SE) and triplet 
even (TE) channels. 



In order to understand the role of the hard core and the tensor interaction, let us 
solve the Schrodinger equation for the deuteron by using the AV8' nucleon-nucleon 
interaction. The wave function of the deuteron is written as 



f^, = u{r)[Yo{r)(E)Xi{(yia2)]iM + w{r)[Y2{f)(g)Xi{cTia2)]iM 



(1) 



The deuteron wave function is written by the 5-wave and (i-wave components. The 
tensor interaction mixes these two components. In Fig. |6] and Table |2] we show the 
deuteron properties for the wave functions and the various energy contributions and 
radii. In the wave function, the i-wave component is dominant and shows the long- 



S-wave — 
D-wave — 




Fig. 6 Deuteron s-wave and d-wave functions 
obtained using the AV8' nucleon-nucleon po- 
tential. 



Table 2 Deuteron properties using 
the AV8' nucleon-nucleon potential. 



Energy 


-2.24 [MeV] 


Kinetic 
(SS) 
(DD) 


19.88 
11.31 
8.57 


Central 
(SS) 
(DD) 


-4.46 
-3.96 
-0.50 


Tensor 
(SD) 
(DD) 


-16.64 
-18.93 
2.29 


LS 


-1.02 


P(Z)) 


5.78 [%] 


Radius 
(SS) 
(DD) 


1.96 [fm] 
2.00 [fm] 
1.22 [fm] 
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tail due to the weak binding of 2 MeV. In the short-range part less than the 0.5 fm 
region, the i-wave function is largely reduced due to the short-range repulsion in 
the central interaction. Looking at the li-wave component, its amplitude starts from 
the origin, because of the centrifugal barrier in the L = 2 partial wave, which is 
absent in the s-wave case. These features of the s-wave and tZ-wave components have 
to be expressed using the shell model framework in finite nuclei as the deuteron- 
like tensor correlation. These features will be treated for finite nuclei by using the 
tensor optimized shell model (TOSM) and the unitary correlation operator method 
(UCOM) to be discussed later 

Among the expectation values of all the two-body interactions in the deuteron, 
the tensor interaction has the largest contribution of about — 17 MeV. This expecta- 
tion value is four times of that of the central interaction. This tensor interaction is 
the origin of the J-wave mixing in the wave function. From these results, it is found 
that the dominant energy contribution comes from the coupling of the s- and d- 
wave components by changing the relative orbital angular momentum by 2, AL ~ 2. 
As for the radius of the deuteron, we decompose the radius into the i-wave and d- 
wave components and normalize them by using the corresponding amplitudes. It is 
interesting to see the size difference between two angular components, where the 
cf-wave size is much smaller than the i-wave size. This compact li-wave structure 
produces high momentum component caused by the tensor interaction, namely the 
pion exchange effect. Hence, we can learn the role and the properties of the tensor 
interaction in finite nuclei from the deuteron. That is, the tensor interaction creates 
the relative d-wa\e component in the wave function, which is spatially compact and 
involves high momentum components. These features should appear in finite nuclei 
as the deuteron-like tensor correlation. We shall treat this correlation in terms of the 
tensor optimized shell model (TOSM) in finite nuclei. 

As for the di-neutron correlation, we have a moderate intermediate attraction with 
a short range repulsion as shown in Fig. |5] There is no tensor interaction and the 
relative motion is completely described by the central interaction. We are aware that 
there is no bound state in the '^o channel, but that the scattering length is negatively 
very large a = —18.5 ±0.4 fm ^. This negatively large scattering length indicates 
that the di-neutron system is close to develop a bound state. Hence, for a system like 
"Li, we expect a strong di-neutron clustering phenomenon in the halo region. For 
the quantitative account we ought to use the A^A^ interaction for this phenomenon. 



1.4 Wave functions for ^'^^Liand ^^Li 

We write the wave functions of the Li isotopes in order to understand the standard 
shell model state, the di-neutron clustering and the deuteron-like tensor correlation. 
It is illustrative to start writing the ^Li wave function. 
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\'Li) = Ci|(.i/2)2 (51/2)2 (p3/2);,(p3/2)t)y=3/2 (2) 

+ C'2|(ii/2);i(il/2)v(P3/2)?i(P3/2)v7=o(/'l/2)vy=0>y=3/2 
+ C3|[('Sl/2);r(il/2)v]7=l(P3/2);i(/'3/2)t[(/'l/2);r(Pl/2)v]y=l)7=3/2 



We have written here only the dominant components explicitly where n and v for 
each configuration denote proton and neutron, respectively. The term with the ampli- 
tude C\ corresponds to the standard shell model state. The term with the amplitude 
C2 corresponds to the main component of the two neutron pairing states, where a 
two-neutron pair couples to J^ =Q^ . The term with the amplitude C3 corresponds to 
the main component of the deuteron-like tensor correlation states, where a proton- 
neutron pair couples to J^ = 1+. 

The di-neutron clustering correlation, which is associated with the C2 amplitude 
component, should involve further particle states in sd and higher shells. As for 
^Li, the di-neutron clustering correlation provides a similar structure as the BCS 
state due to the fact that the nuclear density of the surface neutrons is ordinary 
as expected from the standard size neutron separation energies listed in Table. 1. 
With the increase of the neutron number, the nuclear density of the surface neutrons 
becomes very small and hence the di-neutron clustering correlation should show up. 
This change of the di-neutron clustering correlation due to the nuclear density is 
related with the BCS-BEC crossover. On the other hand, the deuteron-like tensor 
correlation, which is associated with the C3 amplitude component, needs excitation 
of a proton-neutron pair with J^ — \^ from occupied states to unoccupied states. 
We have to include each particle state of the proton-neutron pair up to very high 
angular momentum state. 

We write the '"Li wave function in terms of the ^Li wave function. 

\'^Li)^^[\'^Li)-K\X„)] (3) 

The additional neutron may enter the pj 72 orbit in the shell model state. The ad- 
dition of one more neutron to the pair correlated state has an effect to weaken the 
pairing correlation due to the blocking effect of the neutron. The addition of one 
more neutron to the deuteron-like configuration is very interesting, since the addi- 
tional neutron may go into the p-^n orbit or into the ij /2 orbit in the shell model 
state. If the last neutron goes into the pi « orbit, the deuteron-like correlation is 
weakened by the additional neutron. Instead, if the last neutron goes into the sin 
orbit, the deuteron-like correlation is not weakened, because 2p-2h states with the 
use of i] rt orbit are not important for the deuteron-like correlation. Hence, there 
should appear the competition of the neutron 51/2 and pii2 configurations in '"Li. 
We write the "Li wave function in terms of ^Li wave function. 

|"L/)=^[|^L/)x|;t™)] (4) 
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In this case, there is an important physics to be added in addition to all the inter- 
esting phenomena in '"Li. As for the pair correlated state, two neutrons block the 
pair correlated state by entering in the pi m orbit. More interesting is the case of 
the two neutrons going into the higher shell orbits. In this case, the two neutrons 
stay in a low density region far from the ^Li core and hence the di-neutron clus- 
tering phenomenon is expected. How large is the attraction due to the di-neutron 
clustering effect for the "Li binding energy needs full account of all the effects. 
As for the C3 component, two neutrons going into the p\j2 orbit generate a strong 
blocking effect of the deuteron-like correlation and this configuration is disfavored 
by the tensor interaction. On the other hand, when the two neutrons go into the S\j2 
orbit, the deuteron-like correlation is not disturbed and therefore this configuration 
is favored. Hence, as the consequence of these di-neutron clustering correlation and 
the deuteron-like tensor correlation, these two correlations cooperate to wash out 
the N=8 magic structure and provide the interesting halo phenomenon in "Li. 



2 Di-neutron clustering and the hybrid-Vr model 

In "Li lfT4l and ^He ifTSl . abnormally large matter radii were observed experimen- 
tally. This phenomenon was interpreted as a result of the halo structure, where two 
valence neutrons are spatially extended around the core nucleus due to their weak 
binding. It is important to investigate the dynamics of the motion of valence neu- 
trons for the understanding of the halo structure. In such a situation, it is necessary 
to develop a theoretical method to handle spatially extended structure of the Bor- 
romean system consisting of core nucleus and two neutrons. For this purpose, we 
have developed the hybrid-VT model as the most suitable model to describe the halo 
structure. In the hybrid-Vr model for two-neutron halo nuclei, shown in Fig. |7] the 
mean field nature of each valence neutron can be described in the V-type basis states 
(the cluster orbital shell model (COSM), Fig.|2](a)). Further, the explicit neutron- 
neutron correlation is treated in the T-type (Fig.|7](b)) basis states. Here, the basic 
idea of this hybrid-Vr model is presented, and in the next sub-section we explain 
the formulation of the hybrid-Vr model in detail. 

In the weak binding system it is necessary to consider the large spatial extension 
of single particle wave functions of valence neutrons. This situation corresponds to 
the coupling of the valence neutrons to continuum states. Suzuki and Ikeda proposed 
the cluster orbital shell model (COSM, V-type) ||T6l, and appUed to ^He and "Li 
|[ni[l8l|T9l|20l|21]l22l. The COSM is one extension of the shell model, in which the 
spatially extended character of valence neutrons can be treated. The V-type coordi- 
nates in the COSM are the suitable coordinates to express the mean field property of 
the valence neutrons, and therefore can express the shell model properties of ^'^He 
and '""Li most effectively. However, from the analysis of neutron-rich nuclei with 
the COSM, it was shown that the binding energies of the Borromean nuclei can 
not be described quantitatively. Furthermore, it was found that in a weakly bound 
system, the n-n clustering correlation, namely the di-neutron clustering correlation. 
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Fig. 7 The coordinates of the hybrid-VT model and their relation with those of the Jacobi coor- 
dinates. The V-type coordinates in (a) and the T-type coordinates in (b) are shown in the upper 
figures and the corresponding ones in the Jacobi coordinates are shown in the lower figures. 



becomes important to provide an extra binding energy. This is characterized by the 
participation of many J'' = 0+ pair configurations with large single particle orbital 
angular momenta 123] |24l . 

Hence, it is important to include the physical effect of di-neutron clustering cor- 
relation explicitly in two-neutron halo nuclei, and the T-type wave function is suit- 
able for this purpose ITSl l22l l23l l25l . We combine the T-type basis functions with 
the V-type ones as shown in Fig. |7] It has been shown that this hybrid-Vr model 
describes the di-neutron clustering correlation with full convergence of the binding 
energy and radius as will be discussed in the following sub-sections. Hence, the 
hybrid-Vr model wave function involves two-kinds of physical correlations as the 
mean-field and di-neutron clustering correlations, in the two-neutron halo nuclei. 
The hybrid-Vr model is also a special case of few-body systems; a three-body sys- 
tem is generally described by using the Jacobi-coordinates ||26l as shown in Fig.|2] 
In the case of coie+n+n, the core nucleus has a large mass in comparison with that 
of valence neutrons. Therefore, as shown in Fig. Ill the V-type coordinates (771,772) 
correspond to the symmetric F-type Jacobi coordinates, {X,x) and {Y,y). 

In this section, we explain the construction of the COSM and the hybrid-Vr 
model ||23|24l|25l|27l|^|23[30l|ll]|32l|33l|33[35l[3aandits appUcation to 
the three-body systems: ^He and '^Li. The subject of the di-neutron clustering is 
discussed in terms of BCS-BEC cross over in recent literature ll37l l38l [39l . 
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2.1 Formulation ofhybrid-VT model 

We consider the hybrid-Vr model for a spatially extended core+n+n system. Here, 
we derive the three-body Hamiltonian from the A-nucleon system. This considera- 
tion is useful when we extend this model to include the core excitation in the halo 
nuclear system. The A-body Hamiltonian is given as 

A A 

H^Vti- T,,n + y V.J = T + V, (5) 



A 


A 


E^<- 


- Tan + Yj ^'ij 


(=1 


i>i 



where T = L^/ ^ Tan is the kinetic energy operator of the system after removing 
the center-of-mass motion {Tm) and V = Y^Vij is the two-body potential energy. 
We decompose the Hamiltonian of an A-nucleon system into a core part with Ac 
nucleons and A^(= A —Ac) valence neutrons. The relative coordinates between the 
core and the valence nucleons are given as rj,- = r, — j-Y.i=\ i"/, as shown in Fig.|7] 
The kinetic energy term T is rewritten as 



A 

I 

i=\ 



r = £f,-7;,„ (6) 



Here, Tc = Y^jli f; — T^,„ is the kinetic energy with r^!j„ being the center of mass 
motion of the core nucleus. The operator p=— //iV^ is the momentum conjugate to 
rj and ji =Ac/{Ac + l)m is the reduced mass between the core and a single neutron. 
The term of p, • p j is the recoil motion from the center of mass system. The potential 
term is similarly decomposed as 

Ac N Ac N 

N N 

= Vc + Y,V, + '£vij, (9) 

!=1 i<j 

where the mean field potential Vj for each valence neutron is given as 

Ac 

V.^J^Vij. (10) 

;=i 

Here, Vc are the potential term of the core nucleus, and the Hamiltonian is rewritten 
as 
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H 



N ^ N 



N N 

(=1 i<j 



Hc + Y. 



N.p2 



2fi 



N r- 



+L 

'<./ 



Vii + 



Pi Pj 



{Ac + 1)^1 



(11) 

(12) 



where the first term H^ = T^ + Vc is the Hamihonian of the core and the second and 
third terms are for valence neutrons. The second term is the single particle Hamilto- 
nian for the relative motion between the single neutron and the core. This defines the 
orbitals for the valence nucleons. The third term is the two-body operator between 
the valence neutrons, which produces the coupling between valence neutrons, such 
as the di-neutron correlation. 

We start with the core part 'P{Ac) and write the Schrodinger equation as 



<PiAc)=YCaMAc), 
a 

//,0(A,) = £,0(A,) , 



(13) 
(14) 



where the index a is the label to distinguish various configurations of the core nu- 
cleus and the amplitudes Ca are those used in Eq. (O. They are determined by the 
variational equations obtained by the energy minimization of Ec- We employ the 
shell model like basis wave function for ^a (Ac)- 

The wave function of the two valence neutrons xi™'') in the hybrid-Vr model 
of two neutron halo nuclei (A^ = 2) is expressed as the superposition of the COSM 
(V-type) and T-type wave functions as 



X{nn) = Xv{^v)+XT{^T) 



(15) 



Here, the coordinate sets ^y and ^j. represent V-type and T-type ones, respectively, 
as shown in Fig. |7] We take antisymmetrization between two neutrons, explicitly. 
The radial components of the relative wave functions are expanded with a finite 
number of Gaussian functions centered at the origin with various length parameters 

ioi. 

In the hybrid-V r model, the total wave function of the A-nucleon system and the 
corresponding Schrodinger equation are given as 



//f(A) = E'ViA) , 



(16) 
(17) 



where the total Hamiltonian H is given in Eq. (fT2l i. We omit the angular momentum 
coupling between the core nucleus and the valence neutrons for simplicity. The op- 
erator £/ is the antisymmetrizer between core nucleons and valence neutrons. The 
mixing amplitudes of the core configurations a are included in the wave functions 
Xa{nn). Eq. ( fTSI l is useful to understand the asymptotic condition of the wave func- 
tion, in which some of the valence neutrons are located far away from the core, such 
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as the tail part of halo structure and the scattering states. This will be discussed later 
in the numerical results of "Li. 

In order to solve Eq. ( fTTI i. we employ the orthogonality condition model (OCM) 
iHllllJlllll instead of the resonating group method (RGM) EH. In the OCM, the 
antisymmetrizer jz/ between core nucleons and valence nucleons is replaced by in- 
troducing the projection operator to remove the Pauli forbidden states from the rela- 
tive motion of the valence neutrons. The projection is expressed by introducing the 
following one-body term in the original Hamiltonian in Eq. (fT2]) . 



,.PF 



iV,- 
k 



''M'\ 



(18) 



where the indices / and k are the labels representing each valence neutron and each 
Pauli-forbidden state for one valence neutron. A^, is the number of Pauli-forbidden 
states for one valence neutron. We take a sufficiently large value for X in the numer- 
ical calculation. 

In the calculation of the matrix element of the Hamiltonian in Eq. ( fT2] i. we fold 
the Hamiltonian by using the wave function of the core nucleus. In the coupled 
channel OCM with the hybrid-Vr model, we obtain the following equation for the 
valence neutrons Xa{nn), 



^ai3+L{55„^+V.^/3+vf%4 



/=1 



|:>'-^+(£^}^«^ 



X^{nn)=EXa{nn), (19) 



where 



Hi 



(0a(A,)|//c#/3(A.)), 



\^^„^ = (0«(A,)|V,#i3(A.)) - (0a(A,)|£v,-,-|</.^(A,)). 



(20) 
(21) 



We explain here the Gaussian expansion method to describe the wave functions 
of valence neutrons Xa{nn) in the hybrid-VT model. The spatial part of the basis 
functions for one relative motion r is given by the following Gaussian wave func- 
tions. 



V/f(r)=A^,(Z7)r'exp(-^)y,(f) 



Ni{b) = 



2/,-{2/+3) 

r{l + 3/2) 



(22) 
(23) 
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where / is the orbital angular momentum. The set of the length parameter b is usu- 
ally chosen in geometric progression ||45| . In V- and T-type basis functions, we 
commonly use these basis states. We expand each relative motion of the hybrid-V T 
model with a finite number of the above basis functions. The Gaussian expansion 
method is able to describe the halo structure very nicely. 

For the coupling with intrinsic spin of neutrons, in the COSM (V-type), we adopt 
the j-j coupling scheme in a sense of the shell model. This representation is suitable 
to express the motion of each valence neutron in the mean field potential provided 
by the core nucleus. In the T-type basis function, we take the L-S coupling scheme. 
This is because the di-neutron pair is considered to have the dominant component of 
the 0+ state with spin singlet state. In the T-type basis, we directly take into account 
the ^Sq component of the two neutrons. Actually, the analysis of ^He provides more 
than 80% in the spin singlet states. In the COSM (V-type basis), its basis wave 
function corresponding the configuration of the core nucleus 0a is given as 



p 

p = {biMM,hJi,h,J} , 



(24) 



where j2/i2 is an antisymmetrizer between two valence neutrons, and C„ y are vari- 
ational coefficients for the basis set. 

The T-type basis function is similarly described as 



jXs 



(25) 



xi.T{nn)=Y^Clj^n V/f ^ (r) , ^^ (R 

q - {brM,l,L,S,J], Zs - [Xxl2^Xxl2]s ■ 

In the calculation of the hybrid-Vr model, we take various sets of orbital angular 
momenta and spins until we reach the convergence of the solutions. The variation 
of the total energy E with respect to the total wave function f (A) is given by 

which leads to the following equations: 

— — — - = ^ — ■ — - = ('27') 

BC'' ' BC'' 



2.2 Application of the hybrid-VT model to ^He 

We discuss the results of the hybrid-Vr model for the two-neutron halo nuclei ^He 
as '^He+n+n and "Li as ^Li+n+n. An extension to a core plus many valence neutrons 
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(e.g. ^He=^He+«+«+«) is straightforward P31l47l . For ^He case, we take a single 
configuration for the ""^He core. The Hamiltonian in Eq. ( fT9] l can be written as 



H = Hr 



m 



K^TTj+vp} 



I ^',77^ + 



Pi P2 

(A,+ 1)M- 



(28) 



Here, He can be replaced by the observed energy of ^He (—28.3 MeV) and we 
can discuss only the relative motion of valence two neutrons in ^He. For the ^He-n 
potential V^,j we adopt the microscopic KKNN potential 1231 |48l . which reproduces 
the observed phase shift between "^He and n. The Minnesota interaction [|49| is used 
for v„n between two valence neutrons, where the exchange mixture u is chosen to be 
0.95. These choices are the same as those in Refs. 1251 . l3Tl and l50l . 

It is important to understand the model performance for ^He as a simple system 
for the description of more complicated systems as "Li. The binding energy of 
the ^He ground state is shown as a function of the channel number in Fig. [8] The 
calculated results in the COSM and the hybrid-Vr model are shown by the open 
and solid circles, respectively. In both cases, the convergence was achieved and 
it is found that the hybrid-Vr model converges much faster than the COSM. In 
the hybrid-Vr model, it is sufficient to take shell model states up to the rfj ;2-shell 
orbit (Channel number is five) using the V-type coordinates and add the di-neutron 
channel of Z = L = states in Eq. ( l25T l using the T-type coordinates. Thus, the di- 
neutron correlation is very important for the ground state of ^He. The hybrid-Vr 
model is an efficient framework to treat the di-neutron correlation in the shell model 
basis. 

The radii of ^He are converged as 2.46 fm and 2.06 fm for matter and charge 
rms radii, respectively. These values agree with the recent observations |[T]|52l. The 
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Fig. 8 The binding energy of the ''He ground state as a function of the channel number 1231 . The 
open circles are the results of the COSM model and the solid circles ai'e those of the hybrid-Vr 
model. 
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core-n l^"' 



Fig. 9 The two neutron correlation density p„„{''core-ni S) in ^He in the r and 9 plane [51|. Here, 
I'core-ii denotes the relative distance between the core and one of the two neutrons and 9 the opening 
angle between the two neutrons. 



core-« and n-n mean distances are obtained as 3.42 fm and 4.90 fm, respectively. 
In addition to the radius, the spatial correlations of the halo neutrons in ^He are 
interesting ifTOl [37l l53l l54l and the calculated results are shown in Fig.|9l We show 
the density distribution of halo neutrons p„„ (r, 0) in "^He as a contour in the plane of 
^He-« distance r and the opening angle between two neutrons with the following 
definition 15111591 ; 

p„„(r,0)= f dr'p„„ir,r\e) (29) 

Jo 

p„„(r,/,0) = Sn^r^r'^sinO (f ■'(*'He,r,r', 0)|5(r- r")5(r' - /") 

X 5{e--e')\V\^Ue,r",r"',e')), (30) 

where the total wave function of ^He has three variables, the ^He-n distances with 
r and r' for each neutron and their opening angle 9 in Eq. ( l30l l. We integrate out 
only the variables of the ket part. It is confirmed that the di-neutron type configura- 
tion (a large r and a small 9) gives a maximum value of the density, although the 
density of neutrons is widely distributed. There is another component of the cigar 
type configuration (a small r and a large 9), which coexists with the di-neutron type 
configuration. The characteristics of these two structures in the density distribution 
come from the p^ /j configuration of two neutrons in ^He, which is the dominant 
component by 90.2 % in the ground state wave function. The mixing of the higher 
orbital configurations makes an spatial extension of the distribution and enhances 
the di-neutron component, such as pj ,2 with 4.3%, l^'j /2 with 1.2%, dhr, with 2.6% 
and dy2 with 0.9%. The results indicate that the j-j coupling scheme is well estab- 
lished in ^He. 
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2.3 Hybrid-VT model on di-neutron clustering in ^^Li 

We consider the three-body problem of ^Li+n+n using the orthogonality condition 
model (OCM). The Hamiltonian consists of the similar form as given in Eq. ( |28] | 
for ^He. The difference from the ^He model is the configuration mixing for the ^Li 
core nucleus, because of the small neutron separation energy of ^Li in comparison 
with '*He. This means that we take into account the core excitation in "Li. In this 
section, we first take into account the neutron 0+ pairing correlation of the ^Li core 
and examine this effect on the structures of "Li and "^Li. Later, we include the 
tensor correlation in ^Li. 

Before showing the numerical results, we generally formulate the coupled ^Li+n+n 
model of ' ' Li, in which the configuration mixing is performed for the ^Li core. This 
framework is straightforward to apply when the tensor correlation is included in the 
core part, later. In the coupled hybrid-Vr model of ^Li+n+n, we consider the Pauli 
forbidden (PF) states in the '^Li-n relative motion ll57l . In this model, PF states re- 
moved from the the relative motion depend on the configuration of ^Li, namely the 
orbits occupied by neutrons in the ^Li core. The main configurations are given in 
Eq. ^ and the amplitudes are written as Ci, C2 and C3. The PF states corresponding 
to those wave functions are given as 

{0*1/2, Op3/2 for Ci 

Osi/2, O/73/2, (0;?i/2)vv for C2 (31) 

0*1/2, O/73/2, (0/?i/2);rv for C3 . 

In case of C2 the PF pi/2 orbit is used by the pairing state and indicated as (;?i/2)vv, 
while in case of C3 the PF pi/2 orbit is used by the deuteron-like tensor correlation 
and indicated as {pi/2)7:v- From Eq. ( fTSI ), the wave function of "Li is given as 

f'^("Li) =f ^{[(/)r^zi(««)]'}- (32) 

a 

Here, Xa{nn) represents the wave functions of two valence neutrons, and j and/ are 
the spin of two valence neutrons and the total spin of "Li, respectively. The three- 
body eigenstates are obtained by solving the eigenvalue problem for the coupled- 
channel Hamiltonian given in Eq. (fT2l) . 

//("Li)f("Li) =£("Li)f("Li) (33) 

We discuss the coupling between the ^Li configurations (^^ and the motion of 
valence neutrons. In "Li, the amplitudes Ca of each configuration of ^Li in Eq. ( flJl l 
are determined variationally. Asymptotically, when the two valence neutrons are far 
away from ^Li, the wave function of "Li becomes 
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Xk{nn) > Ca-X'{nn), (34) 



^•'(''Li) 



171. '72- _ 

^3/2 



'71-172- 



J^CarJ' ],X'{nn) 



(35) 



The first equation implies that the asymptotic wave function of the two valence 
neutrons is decomposed into the internal amplitude Ca of the ^Li configuration and 
the relative wave function x^{nn), which is independent of the ^Li configuration. 
This means that the coupling between the valence neutrons and ^Li disappears. As 
for the di-neutron wave function X^(nn), the correlation between the two neutrons 
disappears also at far distance, because the two neutrons do not form bound state in 
the free space. The mixing amplitudes {Ca} of ^Li in Eq. (|35] | are the same as those 
of the isolated ^Li in Eq. ( fT3] ). Contrastingly, when the two valence neutrons are 
close to the '^Li core, the motions of the two valence neutrons dynamically couple 
to the configuration of '^Li in order to satisfy the Pauli principle, which changes the 
mixing amplitudes {Co} in ^Li from those of the isolated ^Li core. 

We now carry out the coupled-channel three-body calculation for "Li. In the 
model, the ^Li-n interaction Vcn is taken as a folding-type potential with the MHN 
interaction P3l [55l l56l . which is constructed from the G-matrix using the bare 
nucleon-nucleon interaction. The folding potential for ^Li-« includes the coupling 
between intrinsic spins of the valence neutron and ^Li (3/2^). This coupling pro- 
duces splittings of the energy levels, for instance 1+-2+ (for the pi ^-neutron) and 
1^-2^ (for the ^i ^-neutron) in the '"Li spectra. 

The important points we wish to study in this calculation are whether the present 
model can solve the under-binding problem and describe the halo structure. This 
is because the three-body model of "Li by using only the inert core model of ^Li, 
which corresponds to the use of the C\ term alone in Eq. (2), does not make a 
bound state ||251 . We also want to see how the pairing correlations act on the binding 
mechanism. The results are shown in Fig.[TO] The binding energy of the "Li ground 
state measured from the three-body threshold is obtained as 0.5 MeV by considering 
the pairing correlation in ^Li. The matter radius is obtained as 2.69 fm, which is 
smaller than the experimental value H] |58l . 

We discuss the role of the pairing correlation between valence neutrons in "Li. In 
Fig. [To] two kinds of the energy convergence of "Li are plotted as functions of the 
channel number of the j''^ = 0+ pairing configuration for valence neutrons. One of 
them is the calculation employing only the COSM basis and the other is that with the 
hybrid-Vr basis. In the calculation, we take the first channel as {pin)^, and the or- 
der in which channels are added to the first one is (*i 72)^, (^3/2)^5 (^^5/2)^' {d^iif', 
■■■ , {IjY- The maximum number of channel is 31, where the orbital angular mo- 
mentum and the spin of one valence neutron are / = 15 and j = ^■■ We see 
rapid convergence of the energy in the hybrid-Vr model. This result indicates that 
the pairing correlation between valence neutrons is important to reproduce the weak 
binding state of "Li. This result is similar to the ^He case shown in Fig. dHJ. We 
comment here that the binding energy comes out to be 0.5 MeV, which is larger 
than the experimental value of 0.3 MeV. At the same time, the s-wave component 
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is less than 10 % as compared with the experimental value ^ 50%. In the present 
analysis, we determine the ^Li-« interaction to reproduce the 1+ state at 0.42 MeV 
and the virtual i-wave state just at the "^Li-n threshold energy in the "'Li spectrum. 
In this case, two neutrons of ^'Li are slightly overbound, which indicate that the 
pairing correlation of ^Li partially solves the problem of '"'^Li. This over-binding 
property together with other problems will be removed consistently by considering 
the deuteron-like tensor correlation, which pushes up the (pi 72)" state energetically 
close to the {sxjif' state, as will be shown later 

We also discuss the di-neutron correlation in '^Li. The spatial correlations of the 
halo neutrons in '^Li are interesting ifTol [37l l53l l54l and they are shown in Fig.fTTI 
We calculate the density distribution of halo neutrons p„„(r, 0) in "Li as functions 
of the ^Li-« distance r and the opening angle Q between two neutrons. In order to 
present the realistic case, which reproduce the large s^ component of halo neutrons 
||3l . we show here the result of calculations with the tensor optimized shell model 
(TOSM), which includes both the pairing and deuteron-like tensor correlations in 
^Li II59I . This TOSM wave function contains the i-wave component in "Li by a 
large amount 47%. We wiU discuss the details of the TOSM in the following section. 
In the TOSM case (a), it is confirmed that the di-neutron clustering configuration 
gives a maximum value of the density, although the density of neutrons is widely 
distributed. Contrastingly, the Inert Core case (b) with a small s^ component of 4%, 
does not show much enhancement of the di-neutron clustering configuration and the 
cigar type configuration coexists with the di-neutron clustering configuration. This 
feature of the case (b) is similar to ^He ll53l . These two results indicate the role 
of the s^ component on the formation of the di-neutron clustering configuration as 
follows: The s^ component in "Li increases the amplitude of the tail region of two 
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Fig. 10 Convergence of the ' ' Li ground state energy with respect to the channel number in the 
COSM (open circles) and in the hybrid-VT model (solid circles) 1241 . The dotted line represents 
energy to which the calculation converges (—0.50 MeV) in the present model setting. 
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Fig. 11 Two neutron correlation density p„„(?'core-ni 6) for "Li 1591 . The case (a) is the calculation 
with the TOSM of Li and (b) is the Inert Core case of Li, respectively. 



neutrons far from ^Li, and these neutrons tend to come close to each other to gain the 
interaction energy between them. As a result, the di-neutron clustering configuration 
is enhanced, although the spatial distribution of neutrons are still wide. The spatial 
distribution of two neutrons also affects the opening angle 9, where the TOSM case 
having large di-neutron component, shows a smaller 9 value (65 deg.) than the Inert 
Core one (73 deg.). 



3 Continuum and resonance states in complex scaling method 

(CSM) 



The halo nuclei have extremely small binding energies and it is very important to 
take into account the continuum and resonance states for quantitative account of the 
halo nuclei. At the same time, neighboring nuclei have odd numbers of neutrons 
and often those nuclei do not have bound states. Hence, in order to obtain precious 
informations on the neutron-core potentials, it is important to describe resonance 
states of the neighboring nuclei. In this section, we would like to develop a powerful 
method of treating continuum and resonance states as if they are bound states. This 
method is called the complex scaling method (CSM). We can also apply the CSM 
for excitation functions of halo nuclei. We would like to emphasize here that there 
has not been any method to treat three-body unbound systems. In particular, there is 
a case where two-body system out of the three-body system may be in the resonance 
state. We can treat these interesting possibilities in the CSM. 
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3.1 Formulation of CSM 

We explain the CSM, which describes resonances and non-resonant continuum 
states of a many-body system. Hereafter, we refer to non-resonant continuum states 
as simply continuum states. In the CSM, we transform every relative coordinates 
{r,} of the system such as core+n+n model, by the operator Uq as 

Ue ■■ r; -^ r,- e''' for i=l,--- ,N , (36) 



where is a scaling angle and A^ the total number of particles in the system. The 
Hamiltonian H is transformed into the complex-scaled Hamiltonian Hg ~ UqHUq , 
and the corresponding complex-scaled Schrodinger equation is given as 

HgWi = £f/, (37) 

^7^g{3/2),ex^7^|^.^,e|)^ (38) 

where X stands for the number of degrees of freedom. The phase factor ef^/^)'^-^ 
is attached here due to the phase freedom of wave function and originates from 
the Jacobian in the integral over the coordinates. In the three-body model of ^'Li, 
X ^2. The eigenstates fg are obtained by solving the eigenvalue problem of Hg in 
Eq. ( |37] |. In the CSM, we obtain all the energy eigenvalues E of bound and unbound 
states on a complex energy plane, governed by the ABC theorem 1601 . In this theo- 
rem, it is proved that the boundary condition of Gamow resonances is transformed 
to the damping behavior at the asymptotic region. The Gamow resonance is a pole 
of S'-matrix and has an complex energy eigenvalue of £ =Er — ir/2, where Er and 
r are the resonance energies measured from the lowest threshold and the decay 
widths, respectively. 

For simple understanding of this theorem, we consider the asymptotic wave func- 
tions for Gamow states in the two-body case. The Gamow states with complex 
wave number kp are described by the outgoing waves e\p{ikpre'^). It is easily un- 
derstood that the bound state wave functions maintain the damping behavior for 
< n/2. The wave functions of resonances, which had divergent behavior origi- 
nally as e'*^«'" = e'('^'-'r-)'- = ^7-'- . g''^-'", behave as 

^ikRi-e'^ ^ gi(Kr-iyr)re'^ _ ^/(-(k-,— jy, )(cose+;sine) 

— „(-K:rSine+7rCos9)r ((K-,.cos8+7rSine)r /oq\ 

This equation shows that the divergent behavior of the resonant wave functions 
is regularized when we take the scaling angle Q to be larger than the angle Or — 
tan^' (^) of the resonance position Kr — ijr- This damping condition enables us to 
use the same theoretical method to obtain many -body resonance states as those used 
for bound states. For a finite value of 0, every Riemann branch cut is commonly 
rotated down by 20. We can identify the resonance poles of complex eigenvalues 
without any ambiguities. 
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Fig. 12 A schematic distri- 
bution of energy eigenvalues 
of the Borromean coie+n+n 
system in the CSM, where 
the origin of energy is chosen 
as the three-body threshold 
energy l43l . 
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In the wave function, the dependence is included in the variational coefficients 
in Eqs. (l24l i and dZSl ) as {C^Vl ^'^'^ i^a rl' respectively. The wave functions are 
expanded with the finite number of the L^ basis functions, so that all the obtained 
states are discretized on the complex energy plane. The stability of the calculated 
matrix elements of resonant and continuum states using the CSM has been shown in 
many works ll27l l43l 1611 . For continuum states, we adopt the discretized represen- 
tation using the L? integrable basis functions. This discretization has been checked 
to reproduce the genuine continuum states by using the CSM ll62l l63l 1641 . 

In the study, we use the CSM not only to search for the resonance positions, 
but also to calculate the strength functions, such as £1 responses. This is related to 
the continuum level density of the scattering states. We have shown that the CSM 
provides us with the accurate continuum level density even if the states are dis- 
cretized. This fact means that the continuous strength function into scattering states 
can be obtained in the CSM, which is performable in the many-body case. So far, 
we have succeeded to apply this characteristics of the CSM to calculate the electric 
responses, Gamow-Teller strengths, nucleon-removal strength and so on. 

In the calculation of the strength function, we need a complete set of the 
coreH-«H-n system including bound, resonant, and continuum states. We express this 
complete set using the complex-scaled eigenstates fg obtained in the coreH-«H-n 
model. We briefly explain the extended completeness relation (ECR) using the CSM 
II27I l6n l65l . When we take a large 9 like in Fig. 12, three-body scattering states 
are decomposed into three categories of discrete three-body resonances, three-body 
continuum states of coreH-nH-«, and two-body continuum states of [coreH-«]ies-i-«. 
Here, the [coreH-«]res-i-n two-body continuum states are obtained on the branch cuts, 
whose origins are resonance positions of the core-nn system, as shown in Fig. 12. 
Using all the unbound states, we introduce the extended three-body completeness 
relation (ECR) of the complex-scaled Hamiltonian Hq as 
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i=ii^v')(^v'i 

= {three-body bound state} 

+ {three-body resonance} 

+ {three-body continuum states of coreH-nH-«} 

+ {two-body continuum states of [coreH-«]res+n} , (40) 

where { 0^ , <Py } are the complex-scaled wave functions and form a set of biorthog- 
onal bases. This relation is an extension of the two-body ECR 1431 |6T1 . Because 
the detailed definition of the biorthogonal bases is written in the previous works 
ll27ll6Tl . we only briefly explain it here. When the wave number ky of <Py is for 
discrete bound and resonance states, the adjoint wave number kv of 0y is defined as 
kv = — fcy, which leads to the relation <Pv = {<Pv)* EIIMIISI- For continuum states, 
the same relation of the bi-orthogonal states of resonances is adopted, because we 
use a discretized representation. In the coreH-«H-« model, the coreH-2n channel is in- 
cluded in the three-body continuum components of the coreH-«H-« system, because 
two neutrons do not have any bound states or physical resonances. 

We explain how to calculate the strength function 5'(£') using the ECR model. 
The strength S{E) is a function of the real energy of the whole system E. We first 
introduce the Green's function ^(E, 7], ?]'), which is used in the derivation of the 
strength ||27]|28l. The coordinates, t] and Tj', represent the set of rt (i — 1 , • • • ,X) in 
Fig. 7. Here, we introduce the complex-scaled Green's function ^^ (£, J? , ??') as 
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In the derivation from Eq. dTTT i to Eq. (|42] |. we insert the ECR of the whole system 
given in Eq. ( |40] |. The total energy E^ corresponds to the eigen wave function 0y . 
The Q dependence of E^ appears only in the continuum spectra. 

The strength function S{E)a for the arbitrary operator Oa, in which a is the 
quantum number for the operator, is defined using the ordinary Green's function as 

Sa{E) = f{%\0i\<Py){$y\0a\%) S{E-Ey) 

jdndn'%{n)o\':^{E,n,i)Oa%{n') 



Im 



(43) 



For simpUcity, we omit the labels of the angular momenta and their z components of 
the wave functions and of the operators. The wave function % is the initial state. We 
also consider the sum rule value of the strength Sa{E) in Eq. (l43l l. which is defined 
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by the integration of Sa {E) over the real energy E. Using the completeness relation 
of the final states of ''He, the sum rule value is given as 



dE Sa{E) =T{%\0U<l>v){1>v\0a\%) 

= {%\OlOam ■ (44) 

Thus, it is also confirmed that the energy integrated value oiSaiE) satisfies the ex- 
pectation value of the operator Oa for the initial state. In the case of the E 1 transition 
of halo nuclei, the sum rule value corresponds to the relative distance between core 
and the center of mass of valence neutrons. When Oa is an annihilation operator, 
for example, ^He into ^He+«, the integrated value satisfies the associated particle 
number of ^He, namely the number of three valence neutrons. 

To calculate the strength function Sa{E) in Eq. ( |43] |. we operate the complex 
scaling on Sa{E), and use the complex-scaled Green's function of Eq. (|42] | as 

Sa{E) = -ilm[|rf77«fT7'[%*(77)]«(o;)«^«(£,77,77')0^fo''(7]') 



f^SaAE), 



Sa,v{E) = Im 

n 



e-eS 



(45) 
(46) 



In Eq. ( |46] |, the strength function is calculated using the matrix elements {<i>y \a^\ ^^ ) 
It is noted that the function Sa.v{E) is independent of 9 1271147116111641 . This is be- 
cause any matrix elements are obtained independently of in the complex scaling 
method, and also because the state v is uniquely classified according to the ECR de- 
fined in Eq. ( i40l i. As a result, the decomposed strength Sa.v (E) is uniquely obtained. 
Thus, the strength SaiE) is calculated as a function of the real energy E of the nu- 
cleus of interest. When we discuss the structures appearing in Sa{E), it is useful to 
decompose Sa{E) into each component Sa,v{E) by using the complete set of the 
final state v of the whole system. We can categorize v using the ECR in Eq. ( i40t . 
Because of this decomposition of unbound states, we can unambiguously investi- 
gate how much each resonant and continuum state of the whole system exhausts the 
strength. This can be performed in the Coulomb breakup strengths of ^He, "Li and 
"Be and many other three-body systems. 



3.2 Three body resonance and continuum states in He 

We discuss the results of the complex scaling method (CSM) for the halo system 
^He. We describe ^He as a "^He+n+n system in the hybrid- VT model li23ll27l . Here, 
we briefly recapitulate the important properties of the CSM. The Hamiltonian of the 
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model is the same as the one in Ref. ||23l except for an introduction of a three-body 
interaction; 

3 2 

H = Y,ti-TG + Y, {Van.i + vD + Vnn + V^,^, (47) 

(=1 (=1 

where f,- and Tq are kinetic energies of each particle and the center-of-mass of the 
three-body system, respectively. The "^He core cluster is assumed to have the (Oi)'*- 
closed configuration with the length parameter bc=\A fm, which reproduces the 
experimental charge radius of ''^He. The Pauli-forbidden PF state is the Qs orbit in 
the relative motion. The two-body interactions Van and v„„ are given by the micro- 
scopic KKNN potential II48I for '^We.-n and the Minnesota potential B9l for n-n, 
respectively. These potentials well reproduce the low-energy scattering data of each 
two-body system. 

A phenomenological three-body ^He-n-n interaction ¥„„„ is introduced to fit the 
binding energy of the ''He ground state. This is introduced to overcome the small 
underbinding (few hundreds keV) of the ^He ground state with a frozen ^He core 
assumption. By taking into account the excitation or the dissociation of the "^He 
core, this underbinding problem in ^He is beheved to be solved Il67ll68l . We include 
effectively the excitation of the '^He core inside ^He by this three-body interaction 
term. We introduce the three-body interaction assuming a single Gaussian function 

vi„=V3.-^"•'+'•^^ (48) 

Va = -0.218MeV, v = {0.1/ bcf fm-^ . (49) 

Using this Hamiltonian, the present hybrid- VT model reproduces the observed ener- 
gies and decay widths of ^^^He, simultaneously 1691 , namely, the threshold energies 
of the particle emissions for the He isotopes. 

The three-body eigenstates are obtained by solving the eigenvalue problem of 
the complex-scaled Hamiltonian. We use 30 Gaussian basis functions for one radial 
component in order to achieve stabilization of the calculated results for the position 
of resonances, distributions of continuum states and their transition matrix elements. 
The maximum range of Gaussian basis functions is about 40 fm. In Fig. [T3] we 
show the eigenvalue distribution for 1^ states of ^He. This result is obtained by 
diagonalization of the complex-scaled Hamiltonian of the '^He+n+n model with = 
35 deg. From Fig. [13] we see that all eigenvalues are obtained along three lines of 
rotated Riemann cuts corresponding to two two-body and one three-body continuum 
channels. There is no 1^ resonance. Therefore, these results indicate that the 1^ 
unbound states above the ^He+n+n threshold are classified into two-body continuum 
states of ^He(3/2^,l/2^)H-« and three-body continuum states of ^HeH-nH-«. 

In Fig. [141 we show the obtained 2+2 resonances and continuum solutions which 
are decomposed into two- and three-body continuum states, similar to the 1 ^ spec- 
tra. The whole energy levels of ^He are displayed in Fig. [15] We can see a good 
agreement between our calculation and experimental data. The present calculation 
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Fig. 13 Energy eigenvalues 
of 1^ states calculated with 
the CSM where 6 is 35 deg. 
1271 . Squares and triangles 
indicate the two-body contin- 
uum states of 'He(3/2^)-l-n 
and 'He(l/2^)-Fn, respec- 
tively. Circles indicate the 
three-body continuum states 
of '^He+n+n. 



Fig. 14 Energy eigenval- 
ues of 2+ states calculated 
with the CSM gTJ, where 
two solid circles are 2j'"2 
resonances and other marks 
indicate the same meanings as 
those in Fig. [13] 
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does not predict any 1 ^ resonance in the low excitation energy region. This resuh is 
consistent with the experimental results 1691 [TOl . 

In Fig. [16] we show the Coulomb breakup cross section with respect to the ex- 
citation energy of ^He. This cross section is calculated using the £1 strength with 
the equivalent photon method, considering the experimental resolution ITtTI . The 
target is Pb and the incident energy of the ^He projectile is 240 MeV/nucleon. It is 
found that there is a low energy enhancement in the strength at around 1 MeV mea- 
sured from the three-body threshold energy. This energy is just above the two-body 
threshold (0.74 MeV) of the ^He(3/2^)H-« system 1231 . and the cross section grad- 
ually decreases with the excitation energy. This structure of the strength indicates 
the sequential breakup process via the ^He(3/2^)H-« channel l27l . We also compare 
the strength with experiments fTTl. The obtained result fairly reproduces the trend 
of the observed cross section, especially in the low excitation energy region below 
E ^2 MeV. The height and position of the low-energy enhancement in the strength 
agree well with the experimental data. 
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Fig. 15 Energy levels of 
*He. Unit of energies and 
resonance widths is MeV 



Fig. 16 Coulomb breakup 
cross section of ''He 1511. 
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In this section, we have discussed the complex scaling method (CSM) in order to 
treat resonance and continuum states in nuclear composite systems. As an example 
of the usefulness of the CSM, we have discussed the case of ^He as a core+«+« 
three-body system. We can interpret the calculated results nicely and can distinguish 
any structure in the continuum spectrum. As for the application of the CSM to the 
Li isotopes, we defer the detailed discussions after the introduction of the deuteron- 
like tensor correlation. The interpretation of the experimental results further need 
the participation of the si/2 orbit in the wave functions of the Li isotopes. For this, 
we have to introduce the deuteron-like tensor correlation, which is the subject of the 
next section. 



4 Deuteron-like tensor correlation and tensor optimized shell 
model (TOSM) 



In this section, we would like to discuss the role of the deuteron-Uke tensor corre- 
lation on the nuclear structure in ^He before applying the developed method to the 
halo nucleus ' ' Li. The tensor interaction is one of the most important ingredients of 
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the bare nucleon-nucleon interaction and plays the central role for the formation of 
finite nuclei and nuclear matter In this lecture note, we have shown explicitly the 
case of the deuteron in Sec. 2, where the tensor interaction plays the central role to 
provide a strong binding energy through coupling of the i-wave component with the 
J-wave component. Although the tensor interaction is known to be important, the 
standard approach of nuclear many -body problems is to obtain the effective interac- 
tion (G-matrix) by solving the Brueckner equation to include the high momentum 
components in the effective G-matrix interaction. Hence, we are used to treat a well 
behaved central effective interaction and have not faced to treat the tensor interac- 
tion explicitly for shell model states. This is the reason why we have not developed 
a method of treating the tensor interaction explicitly in nuclear physics. 

Recently it became possible to calculate nuclei up to mass around A ^ 12 11721 
|73l [131 using the realistic nucleon-nucleon interaction. The method used for the 
calculation is the Green's function Monte-Carlo method (GFMC) with the use of 
relative nucleon coordinates. This method introduces various correlation functions 
with many variational parameters in the nuclear wave function. In the GFMC, the 
nuclear structures and binding energies were successfully reproduced by including 
the three-body interaction. One big surprise is extremely a large contribution of the 
one pion exchange interaction, which is about 70 ~' 80 % of the whole nucleon- 
nucleon interaction. In principle, they can extend this method to calculate heavier 
nuclei. As for the tensor interaction, it contributes about 50 % of the whole two- 
body matrix element. It is however extremely time consuming even with the present 
computer power. Hence, it is strongly desired to develop a new method to calculate 
nuclei with large nucleon numbers by using the nucleon-nucleon interaction. 

The nucleon-nucleon interaction has distinctive features that there exists the 
strong tensor interaction at intermediate distance caused by pion exchange and 
strong short range repulsive interaction at short distance caused by quark dynam- 
ics. The explicit form of the nucleon-nucleon interaction is presented in Sect. 2. 
Although these two interactions have totally different characters, it is customary to 
adopt the Brueckner Hartree-Fock theory to integrate out the high momentum com- 
ponents on the same footing and use the resulting G-matrix as an effective interac- 
tion in the shell model. In this way, we lose information of the tensor correlation and 
the short range correlation in the shell model wave function. Hence, we search for 
a powerful method to treat explicitly both the tensor interaction and the short range 
interaction for the study of not only light nuclei but also medium and heavy nuclei. 

There have been two important developments for this purpose. One is to find 
out that the tensor interaction is of intermediate range and hence we are able to ex- 
press the tensor correlation in a reasonable shell model space 17411751 . We name 
this method as the Tensor Optimized Shell Model (TOSM), where the nuclear wave 
function is written in terms of the standard shell model state and enough amount of 
two-particle two-hole {2p2h) states. This TOSM formalism is based on the success 
of the parity and charge projection in the treatment of the pion exchange interaction 
11761 [TTI . We have shown that the tensor interaction could be treated properly by 
taking a reasonable amount of multipoles (/ < 5) in the 2p2h wave functions with 
the optimization of the radial parts of the particle states. The other is the Unitary 



30 Kiyomi Ikeda, Takayuki Myo, Kiyoshi Kato and Hiroshi Toki 

Correlation Operator Method (UCOM) for the treatment of the short range corre- 
lation ||78l |79l . The short range repulsive interaction is of very short range and it 
is suited to treat the short range correlation using a unitary transformation and take 
the approximation to use only up to the two-body operators. This approximation is 
justified because the volume associated with the short range correlation is extremely 
small, where more than three nucleons rarely enter into the small volume. This is 
not the case for the tensor correlation, since the tensor interaction is of intermediate 
and long range as discussed by Neff and Feldmeier 1791 . 

Our idea is to combine these two methods, TOSM and UCOM, for our purpose to 
develop a theoretical framework to describe furthermore medium and heavy nuclei 
beyond the light nuclei using the realistic nucleon-nucleon interaction. We can use 
the TOSM for the strong tensor interaction utilizing the intermediate nature caused 
by finite angular momentum of the relative wave function and the UCOM for the 
strong short range interaction utilizing the short range nature. We use completely 
different methods for these two distinctive characters of the nucleon-nucleon inter- 
action. After demonstrating its power we hope to apply a newly developed method, 
which we name TOSMU, to many nuclei. Using the TOSMU, we aim at under- 
standing the roles of the tensor and short range correlations in nuclei using the bare 
nucleon-nucleon interaction. As a good start, we would like to apply the TOSMU 
to ^He. Hence, there are two purposes for this study. One is to see how this method 
works for the treatment of the bare nucleon-nucleon interaction. The other is to 
compare with rigorous calculations to check the accuracy of the results obtained in 
the TOSMU. From this comparison, we can see how far we can describe the short 
range and tensor correlations and to find what we are supposed to do for further 
improvement of the TOSMU in order to solve the nucleus as precisely as possible. 

On the other hand, we would like to develop a theoretical framework to describe 
wave functions in terms of single particle coordinates, which we call a V-coordinate 
method. This V-coordinate method has ability to describe nuclei with many nucleons 
relatively easier than the T-coordinate method. Furthermore, we are able to describe 
the wave function based on the shell model picture and hence it is easier to interpret 
the calculated results in the shell model sense. The difficulty, on the other hand, 
is to express the correlations of the relative motion between two nucleons, which 
are caused by the short range repulsive interaction and the tensor interaction in the 
nucleon-nucleon interaction. We overcome this problem by developing the TOSMU 
to describe the short range and the tensor correlations simultaneously. 



4.1 Formulation of TOSM 

We explain the formulation of the Tensor optimized shell model (TOSM). We shall 
begin with the many -body Hamiltonian with mass number A. 

A A 

H^Y.Ti-Tcm + Y.^ij (50) 
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with 

yu = ^^j+^Jj+^u+^j""- (51) 

Here, T) is the kinetic energy of all the nucleons with r„„ being the center of mass 
kinetic energy. We take the bare nucleon-nucleon interaction for Vij such as the 
AV8' consisting of central {vfX tensor (vfj) and spin-orbit (vf?) terms. The vf-'"'' 
is the Coulomb term. We describe the many-body system with many-body wave 
function, f , by solving the Schrodinger equation //f = E^. In the TOSM, we take 
the V-coordinates to express f . 

We have discussed the property of the deuteron, whose wave function is ex- 
pressed in terms of the s-wave and flf-wave components. We ought to express the 
deuteron-like tensor correlation in finite nuclei. We show the deuteron- hke structure 
is expressed in terms of 2p-2h configurations in the shell model framework by tak- 
ing the case of ^He. The two nucleons in the i-orbit have two components. One is 
the ^Si pair and the other is the ^^o pair The ^^i pair can be expressed as 

|[(05i/2)']im) ^fi=oWr/=oW[^o(?)®Zi(cTiCJ2)]iM ■ (52) 

Hence, the relative wave function is in the i-state. On the other hand, when two 
nucleons are in the p-orbit, we can write the wave function of the relative motion 
being in the triplet even channel as 

\[{0pi/2?hM) - fL=0Wr/=2W[i2(?)®Zl(cTlCJ2)]lM • (53) 

Hence, the relative wave function is in the cZ-state. As discussed for the deuteron 
wave function being written in terms of the s and d wave components, we are able 
to express the deuteron-like wave function in the shell model framework by taking 
both the above two components. Hence, we have to introduce 2p-2h wave functions 
to express the deuteron-like tensor correlation in the shell model basis. There is the 
other p-wave component, {pj,/2)^- In this case, the spin and the angular momentum 
have a stretched configuration and this state is not mixed by the tensor interaction. 
This consideration naturally leads us to introduce the tensor optimized shell model 
(TOSM). It remains for us to check if we can express the deuteron-like tensor cor- 
relation within reasonable amount of multipoles in the TOSM 1741 . As for the short 
range correlation appearing in the i-wave component in the deuteron wave function, 
we ought to take a clever method. For this problem, Feldmeier et al. have demon- 
strated that the unitary correlation operator method (UCOM) can be used to treat 
the short-range correlation 117811791 . 

In the TOSM, the total wave function f is written in terms of a linear combina- 
tion of Op-Oh and 2p-2h wave functions. 

1' = Co\0)+Y,Cp\2p2h}p. (54) 

p 
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Here, the wave function |0) is a shell model wave function and |(0i)^) for ^He. 
\2p2h) represents a 2p2h state with various ranges for the radial wave functions of 
particle states. We can write \2p2h) as 

\2p2h)„ = I [[<(xi)v4^(x2)f ®[v/2^(xi)«?>S:(x2)]']'')a . (55) 

The suffix A of the wave function indicates anti-symmetrization of the wave func- 
tions. Here, p denotes a set of quantum numbers of 2p2h states, which are expressed 
with particle (hole) wave functions i//^ iw'af- The index n is to distinguish the dif- 
ferent radial components of the single-particle wave function i/a. The index a is a 
set of three quantum numbers, /, j and t-, to distinguish the single-particle orbits, 
where I and j are the orbital and total angular momenta of the single-particle states, 
respectively, and t- is the projection of the nucleon isospin. The normalization fac- 
tors of the two particle states are included in the wave functions given in Eq. (|55] |. 
For ^He, the coupled spin, J, of two nucleons is 7 — or 7 = 1 . We omit writing 
the coupled isospin, which should be either or 1 depending on the value of J. We 
have used Gaussian functions for radial wave functions to express more effectively 
compressed radial wave functions 1,741 . The shell model technique is used to calcu- 
late all the necessary matrix elements, which are expressed explicitly in Ref. |80|. 
In more heavier nuclei, such as '^C and '^O, the dominant configuration |0) can be 
extended to the superposed ones, which includes the few h(0 configurations such as 
2Ti(0. This part describes the low-momentum component of the wave function and 
the tensor force contribution is not so decisive to determine the nuclear structure. 
The high momentum part is treated by considering the 2p2h excitations \2p2h)p 
from the each low-momentum configurations. In that case, the spatial shrinkage of 
particle states becomes important. 

We explain the Gaussian expansion technique for single -particle orbits B3ll45l . 
Each Gaussian basis function has the form of a nodeless harmonic oscillator wave 
function (HOWE), except for the \s orbit. When we superpose a sufficient number of 
Gaussian bases with appropriate length parameters, we can fully optimize the radial 
component of every orbit of every configuration with respect to the total Hamilto- 
nian in Eq. ( fSOl l. We construct the following ortho-normalized single-particle wave 
function i//^ with a linear combination of Gaussian bases {^a] with length parame- 
ters ba,m- 

Na 
V/^W = E '^a,'n ^a{r,ba,n) for n = !,••• ,A?„ . (56) 

m=l 

Here, Na is the number of basis functions for a, and m is an index that distinguishes 
the bases with different values of ba,m- The explicit form of the Gaussian basis 
function is expressed as 
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M>:ba,,n) = Niiba,,n} r' e-^''/^-'-''/^ [Y,ir),x^/2]jXu , (57) 

1 



Ni{ba,„) = 



.-(21+3) - ' 
'a.m 



r(/ + 3/2) 



(58) 



The coefficients {£/„,„} are determined by solving the eigenvalue problem for the 
norm matrix of the non orthogonal Gaussian basis set in Eq. ( fSTJ i with the dimension 
Na- Following this procedure, we obtain new single-particle wave functions {vSI 
using Eq. ( |56l ). 

We choose the Gaussian bases for the particle states to be orthogonal to the oc- 
cupied single-particle states, which is O^i « in the ^He case. For Os^n states, we 
employ one Gaussian basis function, namely, the HOWE with the length parame- 
ter bf)s^p,m=\ = bos- For Isin states, we introduce an extended Is basis function 
orthogonal to the Osi/2 states and possessing a length parameter bis^,„ that differs 
from b()s 1741. In the extended Is basis functions, we change the polynomial part 
from the usual Is basis states to satisfy the conditions of the normalization and the 
orthogonaUty to the Qs state. 

Two-body matrix elements in the Hamiltonian are analytically calculated using 
the Gaussian bases B3l l45l . whose explicit forms are given in Ref. fl801 . In the 
numerical calculation of following, we prepare 9 Gaussian functions at most with 
parameters of various ranges to obtain a convergence of the energy. Furthermore, we 
have to take care of the center-of-mass excitations. For this purpose, we use the well- 
tested method of introducing a center-of-mass term in the many-body Hamiltonian 

El [111. 

^cm = A (^ + Iamco^rI^ - lh(o) , (59) 

\2Aot 2 2 / 

A 2 A ^ 

Pern = ^ P; , Rem = 7 X! r/ , « = —72- • (60) 

Here, m and A are the nucleon mass and the mass number, respectively, and bo^ is 
the length parameter of the HOWE for the hole 0.? state. We take a sufficiently large 
coefficient. A, to project out only the lowest HO state for the center-of-mass motion. 
In the numerical calculation, the excitation of the spurious center-of-mass motion is 
suppressed to be less than 10 keV. 

The variation of the energy expectation value with respect to the total wave func- 
tion f (^He) is given by 

'-WW-'^ (61) 



which leads to the following equations; 



d{^\H-E\'¥) d{W\H-E\^) 

Ti = ^ ' ^n = ^ • (62) 




34 Kiyomi Ikeda, Takayuki Myo, Kiyoshi Kato and Hiroshi Toki 

Here, £ is a Lagrange multiplier corresponding to the total energy. The parameters 
{ba,m} for the Gaussian bases appear in non linear forms in the energy expectation 
value. We solve two types of variational equations in the following steps. First, 
fixing all the length parameters ba,m, we solve the linear equation for {Cp} as an 
eigenvalue problem for H with partial waves up to Lmax- We thereby obtain the 
eigenvalue E, which is a function of {/?«,,„}. Next, we try to search various sets of 
the length parameters {ba.m} to find the solution that minimizes the total energy. In 
this wave function, we can describe the spatial shrinkage with an appropriate radial 
form, which is important for the tensor correlation 1741 . 



4.2 Formulation of UCOM 

We employ the UCOM for the short-range correlation. Feldmeier et al. worked out 
a unitary correlation operator in the form Il78ll79l . 

= U^j (63) 

'■<;■ 

with Cij = exp(— / gij). Here, / and j are the indices to distinguish particles. Here, the 
two-body operator gij is a Hermite operator, and hence C is a unitary operator. We 
express the full wave function f in terms of less sophisticated wave function as 
f = C^. Hence, the Schrodinger equation, //f = fif becomes H^ = E^, where 
H = C^HC. If we choose properly the unitary correlator C we are able to solve 
more easily the Schrodinger equation. Moreover, once we obtain 0, we can then 
obtain the full wave function, f , by the unitary transformation f = C^. Since C is 
expressed with a two-body operator in the exponential, it is a many-body operator. 
In the case of the short-range correlation, we are able to truncate modified operators 
at the level of two-body operators 1781 . 

In the actual calculation of the UCOM, we define the operator gjj as 

gij = 2 [Pr-iArtj) + s{rij)pr.ij } , (64) 

where the momentum pr.ij is the radial component of the relative momentum, which 
is conjugate to the relative coordinate r,y. s{rij) is the amount of the shift of the 
relative wave function at the relative coordinate r,j for each nucleon pair. Hereafter, 
we omit the indices / and j for simplicity. We also introduce R+{r) as 

«+('-) dt 

7jT = l, (65) 

which leads to the following relation, 
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dR+{r) s{R+{r)) 



dr s(r) 



(66) 



In the UCOM, we use R+{r) instead of s{r) to use the UCOM prescription. R+{r) 
represents the correlation function to reduce the ampHtude of the short-range part 
of the relative wave function in nuclei and can be determined for four spin-isospin 
channels independently. The expUcit form of the transformation of the operator for 
the relative motion is given as 

. X t 1 1 

c rc^R+{r), c' p^c ^ -j===pr^===, clc = 1, (67) 

c^sc = s, c'Si2C ~ 512, c'v{r)c ~ v(7?+(r)), (68) 

where the operators 1, s and S\2 are the relative orbital angular momentum operator, 
the intrinsic spin operator and the tensor operator, respectively. v(r) is the arbitrary 
function depending on r, such as potential. 

In the calculation using the UCOM, we parametrize ^+(r) in the same manner as 
proposed by Neff-Feldmeier and Roth et al. |f78l l79l [83l . We assume the following 
forms for even and odd channels, respectively. 



n^W = '"+« (^^j exp[-exp(r/i3)], (69) 

R°^^{r) = r+a(l -exp(-r/7))exp[-exp(r/)3)] . (70) 

Here, a, j5, y are the variational parameters to optimize the function R+ (r) and min- 
imize the energy of the system. They are independently determined for four chan- 
nels of the spin-isospin pair In the actual procedure of the variation, once we fix 
the parameters included in R+{r), we solve the eigenvalue problem of the Hamilto- 
nian using Eq. (|62] | and determine the configuration mixing of the shell model-type 
bases. Next, we try to search various sets of the R+{r) parameters to minimize the 
obtained energy. 

In the present framework of the UCOM, we introduce the UCOM function R+{r) 
for each spin-isospin channel and ignore the partial wave dependence of 7?+(r). It is 
generally possible to introduce the partial wave dependence in the UCOM and then 
R+ (r) functions are determined in each relative partial wave in the two-body matrix 
elements. Here, we consider the specific case of this extension of the UCOM by 
taking care of the characteristics of the short-range correlation. One of the simplest 
cases of this extension is the UCOM for only the s-wave relative motion, since all 
the other partial waves / except for i-wave (/ = 0) have r' behavior near the origin, 
where the short-range hard core is extremely large. Hence, this r' behavior largely 
cuts down the effect of the short-range hard core. However, only the i-wave function 
is finite at the origin, and the behavior in the origin is determined by the hard core 
dynamics. In fact, the method used by Feldmeier et al. 1781 is to determine the 
unitary operator to reproduce the short-range behavior of the s-wave relative wave 
function. 
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When we incorporate the 5-wave UCOM (5-UCOM, hereafter) into the TOSM, 
we extract the relative i-wave component in all the two-body matrix elements in the 
TOSM using the V-type basis expanded by the Gaussian functions. For numerical 
calculations, we prepare the completeness relation consisting of the T-type basis 
functions \T) as 

1 =^|7;.)(7;-|, I?;-) = \M¥1]l',Xs]jXt), (71) 

where the T-type basis is expanded by the two coordinates of the relative part r and 
the center of mass part R of two nucleons, which are the set of Jacobi coordinates. 
The orbital angular momenta of each coordinate, r and R, are I and L, respectively. 
It is easy to prepare the i-wave relative part by considering / as zero in the T-type 
basis. We construct the above completeness relation of the T-type basis states by 
diagonalizing the norm matrix expanded by the finite number of Gaussian basis 
functions for two coordinates. In the actual calculation, we use 12 bases for each 
coordinate, with which convergence is achieved. 

We calculate the matrix elements of the arbitrary two-body operator O including 
the S'-UCOM correlator Q using the V-type basis with indices a and j3 . Here, we 
insert the above T-type completeness relation in Eq. ( fTTT i as 

{Va\CldQ\Vp) - Y.'yam ■ {Ti\CldQ\Tj) ■ {Tj\Vp). (72) 

'■;■ 

The matrix element using the T-type base, {Ti\C^ OCs\Tj) , is calculated for the two- 
body kinetic part and the central and tensor interactions. For the kinetic part and the 
central interaction, the matrix elements conserve the relative angular momentum, 
and then we can easily calculate the matrix elements of the transformed operator 
C^OCs- For the tensor interaction, the sd coupling matrix elements are properly 
treated, in which Q is operated on only the s-wave part of the relative motion. In 
this case, the operator Q acts on the i-wave relative Gaussian basis function 0/=o('"), 
which is transformed as 



Cs ^i=oir) = —f^ \fRU7) ^,=oiR-{r)), (73) 

where R-{r) is the inverse transformation of R+{r), namely, R^{R+{r)) = r. The 
matrix elements of the T-type basis function are calculated using the above trans- 
formed wave function. 



4.3 Numerical results of TOSM for ^He 

It is important to understand the origin of the large binding energy of ^He for the 
study of the Li isotopes. Particularly, it is important to develop a method to describe 
the source of the large binding energy in the shell model language. Hence, we de- 
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Table 3 Optimized parameters in R+{r) in TOSM+UCOM for four channels in fm in the present 
work. 



a p 



singlet even 
triplet even 
singlet odd 
triplet odd 



1.32 0.88 0.36 

1.33 0.93 0.41 
1.57 1.26 0.73 
1.18 1.39 0.53 



scribe in detail the structure of '^He in the TOSM and also the role of the UCOM. 
First of all, we determine the UCOM functions R+{r) for the calculation of the 
TOSMU. In the UCOM, we optimize the R+ (r) function by changing the three pa- 
rameters of a, j3 and 7 to search for the energy minimum in the TOSMU. In Table 
[3] the optimized three parameters in the S'-UCOM are listed. In Fig. [17] R+{r) func- 
tions used in the present study are plotted in comparison with the case in Ref. ll83ll . 
For the odd channel, in accordance with the discussion in Refs. 1791 and ||83]| . we 
cannot find the optimum value of R+ (r) in the two-body cluster approximation of 
the UCOM transformation for the Hamiltonian. Hence, we decide to fix the range of 
R+{r), namely, (5 as the same one adopted in Ref. Il83l and optimize a and 7, while 
the variation of R+ (r) for the odd channel does not have significant effects on the 
energy and other properties of '*He in comparison with the original case 1791 [83JI . 
Essentially, two types of parameter set of /?+(r) in the present study and Ref. l83l 
give the similar form of R+ (r) for even channels, in which we omitted the correla- 
tion function for the even channels except for i-waves. This result indicates that the 
correlation functions for the short-range repulsion are uniquely determined for each 
channel. 

Next, we show the calculated results of the energy of "^He as a function of Lrnax 
in Fig. [18] We shall then compare the obtained results with the benchmark calcu- 
lation given in Ref. [[84|. To start with, we show the ordinary UCOM case where 
the UCOM is used for all the partial waves. The calculated results of the energy are 
indicated in Fig. [18] by circles as a function of the maximum angular momentum 
^max- The results show good convergence to reach —19 MeV, while the exact value 
of the few-body calculations is approximately —26 MeV as indicated in Fig. [18] 
Although the binding energy is small, we point out here that we can calculate the 
binding energy directly using the nucleon-nucleon interaction in the TOSMU. The 
tensor interaction matrix element is approximately —50 MeV. On the other hand, in 
the previous study ll74l , we obtained approximately —60 MeV for the tensor interac- 
tion matrix element to check the validity of the TOSM, when we used the G-matrix 
for the central interaction to renormalize the short-range repulsion and retained the 
bare tensor interaction of AV8' in our previous calculation. This fact indicates that 
the treatment of the short-range repulsive interaction is interfering with the contribu- 
tion of the tensor interaction. This is due to a large removal of the short-range part of 
the relative wave functions in the UCOM, in particular, in the (i-wave part of the sd 
coupling of the tensor interaction matrix element, where the tensor interaction pos- 
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sesses some amount of strength. We have also calculated the contributions beyond 
the 2p2h configurations in the TOSM such as 3p3h and 4-p4-h configurations. When 
we include the 4-p4-li configurations within the /9-shell, their contribution to the bind- 
ing energy is approximately 50 keV. This fact denotes that these more complicated 
wave functions contribute very little in the total '^He wave function. 
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Fig. 17 The short range correlation functions, S+(r), for the UCOM in even and odd channels. 
The thin curves are for the R+{r) function of Roth et al. and the thick curves are those of Myo et 
al. 



We have decided to restrict the use of the UCOM to the relative i-wave only for 
the even channel (5-UCOM), where the treatment of the short-range repulsion is 
absolutely necessary. In other partial waves, we have the centrifugal potential that 
cuts out the short-range part from the wave functions of the higher partial waves. 
In this case, we can use the modified interaction and the kinetic energy only for 
the relative i-wave component in the even channels. Since the use of the UCOM 
for the odd partial wave is slightly better, we use the UCOM for all odd partial 
waves. As a starting calculation, we have neglected the 5-UCOM correlation in 
the calculation of the tensor interaction matrix elements. The energy converges to 
—24 MeV, which is now very close to the exact one as shown in Fig. [18] In this 
case, the tensor interaction matrix element is —61 MeV, which becomes close to 
the exact value of —68 MeV. This improvement mainly comes from the increase in 
the sd coupling of the tensor interaction matrix element, however, this calculation 
is still not yet perfect. We have to treat the effects of the short-range repulsion on 
the tensor interaction matrix element. Hence, we have worked out the formulation 
to treat the rigorous s-wave function with the effect of the short-range repulsion for 
the calculation of the tensor interaction matrix elements as explained in the previous 
section. 

The numerical calculation is quite involved in the S'-UCOM case, since the s- 
wave relative wave function with the effect of short-range repulsion should be used 
for the tensor interaction matrix element. We show the calculated results for the 
total energy by the squares in Fig. [18] We see quite a satisfactory result for the 
total energy, which is approximately —22 MeV. We show now all the components 
of the energy for '*He in Fig. [19] All the energy components show the saturation 
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Fig. 18 The energy of '^He in the TOSMU as a function of the the maximum angular momentum 
Anax- The circles are the results where the UCOM is used for all the partial waves. The squares are 
the results using 5-UCOM in the tensor interaction matrix elements. 



behavior as function of Lmax- In the tensor component, the saturation is obtained at 
around Lmax being 8. For the other components, their saturation points are seen at the 
similar Lmax- A very interesting feature is the kinetic energy, which goes up to a large 
value as the tensor interaction matrix element becomes large. As for the comparison 
with the rigorous calculation, we see that Vc satisfies the rigorous value, which is 
approximately —55 MeV. On the other hand, the tensor interaction matrix element, 
Vt converges to —55 MeV, while the rigorous one is —68 MeV. The kinetic energy is 
approximately 90 MeV, while the rigorous one is 102 MeV. The LS matrix element 
is also smaller than the rigorous value. As the net value, the total energy, E, is —22 
MeV and the rigorous value is —26 MeV. A detailed comparison is performed in 
TablelH in which the converged energies in the TOSMU are shown with the rigorous 
calculations. One of the possibilities for the lack of the energy in the TOSMU is 
due to the separate treatment of the short-range and tensor correlations. Although 
the dominant part of the tensor interaction is of intermediate and long range, there 
may remain some small strength in the short-range part of the tensor interaction, 
which can couple with the short-range correlations. This effect can be included by 
extending the truncation of the UCOM transformation in the Hamiltonian to more 
than the two-body level. Three-body term of the UCOM transformation is one of 
the possibilities to overcome the lack of energy in the TOSMU ||78]| . 

In Table |5] we list the mixing probabilities of the dominant configurations in 
^He. The subscripts 00 and 10 represent J and L, the spin and isospin quantum 
numbers, respectively. It is found that the 2p2li configurations with {J, L)=(1,0) for 
the particle pair state are significantly mixed. These spin and isospin are the same 
as those for the deuteron, and thus, this two-nucleon coupling can be understood as 
a deuteron-like correlation 1741 . 
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Fig. 19 Matrix elements of the central interaction (Vc), tensor interaction (Vt) and the spin-orbit 
interaction (Vis) together with the kinetic energy (Kinetic) and total energy (Energy) in the Hamil- 
tonian for ^He as function of L^ax- We observe good convergence for all the matrix elements. 
These values are compared with the benchmark results of Ref. |f84:|, which are indicated by the 
thick short solid lines on the right-hand side of the figure. 



Table 4 Total energy, matrix elements of the Hainiltonian and radius of '^He are compared with the 
benchmark results denoted by FY. Units are in MeV for the total energy and the matrix elements, 
and fm for the radius of *He. 



Energy Kinetic Central Tensor LS Radius 



Present(UCOM) 

Present(5-UCOM) 

FY 



-19.46 88.64 
-22.30 90.50 
-25.94 102.39 



-56.81 -50.05 -1.24 1.555 
-55.71 -54.55 -2.53 1.546 
-55.26 -68.35 -4.72 1.485 



We have been describing '^He as an example of the TOSMU for the strong tensor 
correlation with the use of the UCOM for the short range correlation. This method 
is able to describe the ^He system almost precisely with the use of the bare nucleon- 
nucleon interaction. We are working out further the small difference present between 
the TOSM and the few body many body methods. We believe now that the difference 
comes from the competition of the short range repulsion and the tensor attraction in 
the very short range part of the relative wave function. However, we believe that the 
TOSM is able to describe the deuteron-like tensor correlation and we shall apply 
this method for the description of the Li isotopes. 
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Table 5 Mixing probabilities in the '*He ground state in %. 



(O^-)^o 


82.48 


(o^)ro'(opi/2)io 


2.54 


(0^)ro'[{ 1^-1/2 )(0''3/2)] 10 


2.34 


(0.)ro^[(Op3/2)(0/5/2)]lO 


1.90 


{0s)^i[{0pi/2){0Pi/2)]m 


1.55 


(0.)ro' [(0^5/2 )(0g7/2)] 10 


0.79 


(0.)r„2 (0^3/2 )2o 


0.44 


remaining part 


7.96 



5 Di-neutron clustering and deuteron-like tensor correlation in 
Li isotopes 

The biggest puzzle from the theory side is the large i-wave component for the halo 
neutrons in "Li. If we interpret this fact in the shell model, the shell gap at A' = 8 
has to disappear. However, the mean field treatment of a central force is not able to 
provide the disappearance of the N — di shell gap. So far, there were many theoretical 
studies for "Li 121[10l[l2l[T9l|24l|25][85l[86l[H21[88l[89l|90l and essentially all the 
theoretical works of "Li had to accept that the 1^1/2 single particle state is brought 
down to the Op 1/2 state without knowing its reason [9,,88J. It is therefore the real 
challenge for theoretician to understand this disappearance of the A^ = 8 shell gap, to 
be called s-p shell gap problem, which is discussed in this section by developing a 
framework to treat the deuteron-like tensor correlation explicitly using the nucleon- 
nucleon interaction. The halo structure of "Li is also related with the li-state and 
the Op-state in ^''Li. Several experiments suggest the dual states of the ,?! ^-state 
coupled to the 3/2^ proton state appears close to the threshold energy of ^LiH-« 
together with the dual states of the pi ;2-state ||9Tll92l . 

The tensor interaction, on the other hand, plays an important role in the nuclear 
structure. For example, the contribution of the tensor interaction in the binding of 
"^He is comparable to that of the central force ll84l l93l . The tensor correlation in- 
duced by the tensor interaction was demonstrated important for the "^He-nn system 
II94I I95I l96l . We treated there the tensor interaction in the shell model basis by 2p- 
2h excitations, and found that the {0si/2)^^{0pi/2)^ excitation of the proton-neutron 
pair has a special importance in describing the tensor correlation in '^He 17611961 . 
This 2p-2h excitation causes the strong Pauli-blocking in the '^He+n system for the 
Pi /2-orbit of the last neutron, which contributes to the p-wave doublet splitting in 
^He 1961 . The same effect of the deuteron-like tensor correlation with the Pauli- 
blocking of the additional two neutrons in the pi^ orbit is expected in "Li. The 
occupation of the two neutrons in the pi^ orbit interferes with the deuteron-like 
tensor correlation, which is used to provide a large binding effect of ^He in ^Li. 

Hence, it is important to study the effect of the deuteron-like tensor correlation 
together with the pairing correlation for the s-p shell gap problem in "Li. This is the 
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purpose of this section. We shall perform the configuration mixing based on the shell 
model framework for ^Li to describe the tensor and pairing correlations explicitly. 
In particular, we pay attention to the special features of the tensor correlation. For 
' ' Li, we shall solve the coupled ^Li+n+n problem which treats both correlations and 
investigate further the Coulomb breakup strength of "Li and other observables to 
see the effect of these correlations. 



5.1 Model of Li isotopes 

We shall begin with the introduction of the model for ^Li, whose Hamiltonian is 
given as 

9 9 

HCU)^Y^ti-tG + Y,Vij. (74) 

i=l i<j 

Here, f,, tc, and v,,- are the kinetic energy of each nucleon, the center-of-mass (cm.) 
term and the two-body A^A^ interaction consisting of central, spin-orbit, tensor and 
Coulomb terms, respectively. The wave function of ^Li(3/2^) is described in the 
tensor-optimized shell model 17411961 . We express ^Li by a multi-configuration, 

^C^i)='£Ci<pf/^\ (75) 

i 

where we consider up to the 2p-2h excitations within the Op shell for <Pj in a 
shell model type wave function, and A^ is the configuration number. Based on the 
previous study of the tensor-optimized shell model Il74ll96l , we adopt the spatially 
modified harmonic oscillator wave function (Gaussian function) as a single particle 
orbit and treat the length parameters ba of every orbit a of Os, Opi/2 ^nd Op^/2 
as variational parameters. This variation is shown to be important to optimize the 
tensor correlation Il74ll76ll77ll96l . 

Following the procedure of the tensor-optimized shell model, we solve the vari- 
ational equation for the Hamiltonian of ^Li and determine {C,} in Eq. (iTSl l and the 
length parameters {ba} of three orbits. The variation of the energy expectation value 
with respect to the total wave function f (^Li) is given by 

^ (-fi-f) = ° ' ^^^^ 

which leads to the following equations: 

d{^\HeLi)-E\^) _^^ d{^\HeLi)-E\^) _^_ ^^^^ 

dba ' dCj 
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Here, E is the total energy of ^Li. The parameters {ba} for the Gaussian bases 
appear in non-linear forms in the total energy E. We solve two kinds of variational 
equations in the following steps. First, fixing all the length parameters ba, we solve 
the linear equation for {C,} as an eigenvalue problem for //(^Li). We thereby obtain 
the eigenvalue E, which is a function of {ba}- Next, we try various sets of the length 
parameters {ba} to find the solution which minimizes the energy of ^Li. In this wave 
function, we can optimize the radial form of single particle orbit appropriately so 
as to describe the spatial shrinkage of the particle state, which is important for the 
tensor correlation 17411761 1771 196| . 

For ' ' Li and '''Li, their Hamiltonians are written in terms of ^Li+«+« and ^Li+n, 
respectively, and are given as 

2 2 

//(iiLi) = //(^Li) + £ r, - r^'' + £ y,„,,+y„,„ (78) 



(79) 



where //(Li), T^, 7^ and T^ are the internal Hamiltonian of Li given by 
Eq. (|74l i, the kinetic energies of each cluster ik — Q for ^Li) and the cm. terms 
of three or two cluster systems, respectively. Vcn.k are the ^Li core-« interaction 
{k= 1,2) and V„„ is the interaction between last two neutrons. The wave functions 
of "Li and '"Li with the spin J and /', respectively, are given as 

f^C'Li) = Y^s^ [[cpf-\xHnn)y] . (80) 

f^V''Li)=f^{[^f/'^z/°(«)]^'}. (81) 

We obtain the coupled differential equations for the neutron wave functions x^° (««) 
and X''o (n ) , where Jo and J'q are the spins of the additional neutron part of " Li and 
'"Li, respectively. To obtain the total wave function f ■'("Li) and f"' ('"Li), we 
actually use the orthogonality condition model (OCM) |[T9l l24l l43l to treat the an- 
tisymmetrization between last neutrons and neutrons in ^Li. In OCM, the neutron 
wave functions % ^^ imposed to be orthogonal to the occupied orbits by neutrons 

in ^Li, which depend on the configuration <P^f in Eq. dTSl l. We obtain the fol- 
lowing coupled Schrodinger equations with OCM for the set of the wave functions 

ixHnn)} for "Li and {z/"(n)} for '"Li, where / = 1, • • • ,A?: 
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Ai^X £ l-^aX-^al, (84) 

aG*,('Li) 

where /i/^C^Li) = {^^J^' \H{^U)\(p]'^-' ). T^^f and 7;^f are the total kinetic energies 
consisting of the relative motions for '*Li and '"Li, respectively. A,- is the projection 
operator to remove the Pauli forbidden states 0a from the relative wave functions 
Il50ll57l . where 0a is the occupied single particle wave function of the orbit a in 
^Li defined in Eq. (l3Tl i. This A; depends on the neutron occupied orbits in the con- 
figuration 0- ' of ^Li and plays an essential role to produce the Pauli-blocking 
in "Li and "^Li, explained later This Pauli blocking term A, reduces the pairing 
and the deuteron-like tensor correlations depending on the occupation of the addi- 
tional neutrons in shell model single particle orbits. In '"Li and "Li, the term is 
particularly effective when the neutron or neutrons occupy the pin single particle 
orbit. The value of /L is taken large as 10^ MeV in the present calculation in order 
to project out the components of the Pauli forbidden states into an unphysical en- 
ergy region. Here, we keep the length parameters {ba] of the single particle wave 
functions as those obtained for ^Li. 

We explain the method of treating the orthogonality condition including the 
particle-hole excitations of ^Li in more detail 12411961 . When the neutron orbit in 
^Li is fully occupied, the orthogonality condition for the last neutrons to this orbit is 
given by A,- in Eqs. ( l82b and ( l83b . When neutron orbits in ^Li are partially occupied, 
such as in the 2p-2li states, the last neutrons can occupy these orbits with particular 
probabilities, which are determined by the fractional parentage coefficients of the 
total wave functions of '" "Li consisting of ^Li and the last neutrons. 

We describe the two neutron wave functions x(™) in Eq. ( |82] ) for "Li precisely 
in a few-body approach of the hybrid-Vr model |fT9ll23ll24l : 

xf''{nn) = xf''i^v) + xf''i^T), (85) 

where £, y and ^ j are V-type and T-type coordinate sets of the three-body system, 
respectively, as shown in Fig. [7] 

Here, we discuss the coupUng between ^Li and the last neutrons, whose details 
were already explained in the pairing -blocking case 1241 l28l l50l l97l . We consider 
the case of "Li. In the present three-body model, the Pauli forbidden states for the 
relative motion provides the Pauli-blocking effect caused by the last two neutrons 
l24l[50l . This blocking depends on the relative distance between ^Li and the two 
neutrons, and change the structure of ^Li inside "Li, which is determined varia- 
tionally to minimize the energy of the "Li ground state. Asymptotically, when the 
last two neutrons are far away from ^Li (^y j- — > °°), the effects of antisymmetriza- 
tion and the interaction between ^Li and two neutrons vanish in Eq. ( |82l i. There- 
fore, any coupling between ^Li and two neutrons disappears and ^Li becomes its 
ground state. Namely, the mixing coefficients {C,} are the same as those obtained 
in Eq. ( l75l l. Contrastingly, when the two neutrons are close to ^Li, the two neutrons 
dynamically couple to the configuration <t>- ' of ^Li satisfying the Pauli principle. 
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This coupling changes {C,} of ^Li from those of the ^Li ground state, and makes 
the tensor and pairing correlations to be different from those in the isolated case. 
For '"Li, the similar coupling scheme is considered. The dynamical effect of the 
coupling arising from the PauU-blocking is explained in the results in detail. 



5.2 Effective interactions 

We explain here the interactions employed in Hamiltonians in Eqs. ( l74b . dTSl l and 
( |79l ). Before explaining the present interactions, we give a brief review of the sit- 
uation of the treatment of the effective interactions for the study of ^"' "Li. As 
was mentioned, most theoretical studies based on the three-body model of "Li em- 
ploy the state-dependent ^Li-n potential where only the s-wave potential is made 
deeper than other partial waves |J9l, while the ^Li core is described as inert. This 
state-dependence in the ^Li-« potential is phenomenologically determined in order 
to satisfy the experimental observations of a large s^ component and a two-neutron- 
separation energy of "Li, and a virtual i-state in "'Li, simultaneously. On the other 
hand, for the nn part, the interaction having a mild short-range repulsion 1251 [86l 
or the density-dependent one are often used ifTOl . However, even in the microscopic 
cluster models using an unique effective A^A^ interaction consisting of the central 
and LS forces |[T2l[85l . the s-p shell gap problem in "Li and '"Li cannot be solved 
simultaneously. From these results, we consider that the usual approach based on 
the effective central and LS interactions may be insufficient to explain the exotic 
structures of '""Li. For this problem, even the so-called ab-initio calculations us- 
ing the realistic A^A^ interactions, such as Green's function Monte Cairo 1981 , do not 
provide good results for "Li. 

In this study, we focus on the tensor correlation, which is newly considered to 
figure out the s-p shell gap problem. To do this, we extend the three-body model of 
"Li to incorporate the tensor correlation fully, in particular, for the '^Li part. In the 
present study, our policy for the study of " Li is to use the experimental informations 
and the corresponding theoretical knowledge for ^Li and '"Li as much as possible. 
Following this policy, we explain our interactions in three terms; v,y of //('^Li) in 
Eq. ( l74l i. core-« Va, and n-n V„„ of the Hamiltonians in Eqs. dTSl ) and ( l79l ). 

For the potential V„„ between the last two neutrons, we take a realistic interaction 
AV8' in Eq. (fTsT l. Our interest is to see the n-n correlation in the two-neutron halo 
structure, and therefore it is necessary to solve two-neutron relative motion without 
any assumption. For this purpose, our model space of two neutrons using the hybrid- 
VT model shown in Eq. ( [85] l has no restriction and wide enough to describe the short 
range correlation under the realistic nuclear interaction AV8'. Therefore, there is no 
parameter in the potential V„„. 

The ^Li-7! potential, ¥„,, in Eqs. (iTSl l and ( |79] l is given by folding an effective 
interaction, the MHN interaction 15611991 , which is obtained by the G-matrix calcu- 
lation and frequently used in the cluster study of light nuclei lT9l l36l l43l l50l l56l . In 
the '^Li+n system, the folding potential for the ^Li density calculated by using H.O. 
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wave function has been discussed to reproduce the proper energies of the '"Li spec- 
tra lfT9l l24l |50l . Furthermore, considering the small one-neutron-separation energy 
of ^Li and a long-range exponential tail of the density, we improve the tail behavior 
of the folding potential to have a Yukawa type form Il24ll28l . Any state-dependence 
is not used in the present ^Li-« potential, such as a deeper potential for the 5-wave. 
This is possible because the Pauli blocking effect of the single particle state is in 
action and the state with the pi/2 orbit is pushed up in energy and becomes close 
to the state with the s^n state ll24l |50l . We will discuss the results on "'Li after 
the discussion on "Li. We introduce one parameter 5, which is the second-range 
strength of the MHN G-matrix interaction in the calculation of the ^Li-« potential 
as shown in Table |9l The paramter 8 is to describe the starting energy dependence 
dominantly coming from the tensor interaction in the G-matrix calculation Il43ll56ll . 
In the present calculation, we chose this 5 parameter to reproduce the two-neutron- 
separation energy of "Li as 0.31 MeV after working out the tensor and the pairing 
correlation effects as explained later. It is found that this folding potential also repro- 
duces the positions of the p-wave resonances in '"Li, just above the ^LiH-« threshold 
energy ||9T1 . as shown in the results. 

Now we discuss the choice of the interaction between nucleons in the ^Li core; 
Vij in H(^Ia), where we use the limited shell model wave functions up to the p shell 
for the ^Li core in Eq. (ITSl l. Since our main interest in this work is to investigate 
the role of the tensor interaction and the Pauli blocking effect on the two-neutron 
halo formation, we describe the tensor correlation in addition to the pairing cor- 
relation in the ^Li core. Along this line, recently we have many interesting works 
Il76l iTTl 1 1001 I101II102I . We have also studied the role of the tensor interaction in 
the shell model framework, and proposed the tensor-optimized shell model Il74ll96l . 
As a reliable effective interaction considered from those studies, in this calculation, 
we use GA proposed by Akaishi 1961 |TOT] [ToU for v,;, in Eqs. ^, dTSj and d79] l. 
This effective interaction GA has a term of the tensor interaction obtained from the 
G-matrix calculation using the AV8' reaUstic potential keeping the large momen- 
tum space 1101111021 . Since we limit the shell model space, we increase the tensor 
strength by 50%. This increase of the tensor strength has been studied in the full 
TOSM to provide a quantitative account of the tensor correlation in the limited shell 
model space. In GA, the obtained ^Li wave function in Eq. ( fTSl l shows smaller mat- 
ter radius than the observed one due to the high momentum component produced by 
the tensor correlation 1176117711961 . Hence, we have to adjust the central force, which 
is done by changing the second range of the central force by reducing the strength 
by 2L5% and increasing the range by 0.185 fm to reproduce the observed binding 
energy and the matter radius of ^Li in the same manner as done for ^He 



5.3 "^Li 

We first show the results of the ^Li properties, which give a dynamical influence 
on the motion of last neutrons above the ^Li core in "'"Li. In Fig.|20l we display 
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Fig. 20 (Color online) Energy surface of 'Li with respect to the length parameters ba of Op orbits 
1751 . The two minima indicated by (a) and (b) in the contour map correspond to the states due to 
the tensor correlation and the paring correlation, respectively. 



the energy surface of ^Li as functions of the length parameters of two Op orbits, 
where bo^ is aheady optimized as 1.45 fm. There are two energy minima, (a) and 
(b), which have almost a common ^Opjo value of 1.7-1.8 fm, and a small (0.85 fm) 
and a large (1.8 fm) bop^.^ values, respectively. The properties of two minima are 
listed in Table|6]with the dominant 2p-2li configurations and their probabilities. It is 
found that the minimum (a) shows a large tensor contribution, while the minimum 
(b) does not. Among the 2p-2h configurations, the largest probabilities are given by 

(0'^)h?(OPi/2)?o f*^'" (^^' similar to the results in Ref. f96l, and (0/53/2)01^(0^1/2)01' 
namely the Op shell pairing correlation for (b). These results indicate that the minima 

(a) and (b) represent the different correlations of the tensor and pairing characters, 
respectively. The spatial properties are also different from each other; the tensor cor- 
relation is optimized with spatially shrunk excited nucleons for (a) and the pairing 
correlation is optimized when two Op orbits make a large spatial overlap for (b). 
In Table |6] we show the results of the superposition of minima (a) and (b), named 
as (c), to obtain a ^Li wave function including the tensor and pairing correlations, 
simultaneously. For (c), the favored two configurations in each minimum (a) and 

(b) are still mixed with the Qp-Oh one, and the property of the tensor correlation is 
kept in (c). The superposed ^Li wave function possesses both the tensor and pairing 
correlations. 



5.4 Pauli-blocking effect in Li 



We discuss here the Pauli-blocking effect in '^Li and '"Li. Considering the proper- 
ties of the configuration mixing of ^Li, we discuss the Pauli-blocking effects in '"Li 
and "Li and their difference as shown in Fig.[2T| For (a) in Fig. [21] the ^Li ground 
state (GS) is strongly mixed, in addition to the Op-Oli state, with the 2p-2h states 
caused by the tensor and pairing correlations. 
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Table 6 Properties of 'Li with configuration mixing. The .states (a) and (b) correspond to the each 
energy minimum shown in Fig.|20] respectively. The states (c) is obtained by superposing (a) and 
(b). 







Present 
(a) (b) 


(c) 


Expt. 


E [MeV] 
{Vt) [MeV] 


-43.8 -37.3 
-22.6 -1.8 


-45.3 
-20.7 


-45.3 


R„, Ifm] 


2.30 2.32 


2.31 


2.32±0.02llJ 


Op-O/j 

(OP3/2)oi"(OPl/2 
(0^1/2)r(?(0Pl/2 


)^1 


91.2 60.1 
0.03 37.1 
8.2 1.8 


82.9 
9.0 

7.2 


— 



Let us add one neutron to ^Li for '"Li. For (b) in Fig. [21] when a last neutron 
occupies the Qpin orbit for the p-state of "'Li, the 2p-2h excitation of the pairing 
correlation in ^Li are Pauli-blocked. The tensor correlation is also blocked partially, 
but not fully by the Pauli principle because the Qpin orbit is not fully occupied by a 
last neutron. Accordingly, the correlation energy of ^Li is partially lost inside '"Li. 
Contrastingly, for (c) in Fig. [21] the Is state of '"Li, the Pauli-blocking does not 
occur and ^Li gains its correlation energy fully by the configuration mixing with 
the 2p-2h excitations. Hence, the energy difference between p and s states of '"Li 
becomes small to explain the inversion phenomenon li24l[50l . 

For "Li, let us add two neutrons to ^Li. The similar blocking effect is expected 
for "Li, whose important properties were given in the previous paper IJTSl . For 

(d) in Fig. [2T] when two neutrons occupy the Opi ^-orbit, the 2p-2h excitations 
of the tensor and pairing correlations in ^Li are Pauli-blocked. In particular, the 
blocking of the tensor correlation in "Li is expected to work stronger than the '"Li 
case, due to the presence of the last two neutrons in the /71/2 orbit. Accordingly, 
the correlation energy of ^Li is lost inside "Li stronger than the '"Li case. For 

(e) in Fig.[2Tl (Is)'^ of two neutrons, the Pauli-blocking does not occur, similar to 
the Is state of '"Li. Hence, the relative energy distance between (Op)^ and (li)^ 
configurations of "Li becomes small to break the magicity in "Li. It is found that 
there is a difference in the blocking effects between "Li and '"Li. It is interesting 
to examine how this difference affects the s-p shell gap problem in these nuclei. For 
the pairing correlation, we already pointed out the different blocking effects between 
'"Li and "Li in the previous study li24l . We further consider the blocking effect in 
the dipole excited states of "Li later, which is also different from the ground state 
case. In the previous paper li75l . we examined that the configuration mixing of the 
id-shell for ^Li would give a small influence on the blocking effect on the (li)^ 
configuration of "Li. 
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Fig. 21 Schematic illustration for the Pauli-blocking in "Li. Details are described in the text. 



5.5 '''Li 



We solve '"Li in a coupled ^Li+« model, in which the ^Li-« folding potential is 
determined to reproduce the two-neutron separation energy 82,1 of "Li as 0.31 
MeV. The resonant states are described using the complex scaling method. On 
this condition, we investigate the spectroscopic properties of '"Li. In Table |7] 
it is shown that using the TOSM for ^Li, the dual p-state resonances are ob- 
tained near the ^LiH-« threshold energy. The dual states come from the coupling 
of [{0p-i/2)ji{0pi/2)v]i+/2+^ while the experimental uncertainty is still remaining 
including the spin assignment Q. The 1+ state is predicted at a lower energy than 
the 2+ state due to the attractive effect of the triplet-even ^E channel in the pn in- 
teraction. 



50 



Kiyomi Ikeda, Takayuki Myo, Kiyoshi Kato and Hiroshi Toki 



Table 7 The resonance energies Ey and the decay 
widths r of the p-wave resonance states ( 1 + and 
2+ states) of "^Li in unit of MeV, measured from 
the ^Li+n threshold. The scattering lengths Us of 
the j-states (1^ and 2^ states) are shown in unit 
of fm. We show here the two kinds of the results 
using TOSM and Inert Core for 'Li. 



(£„r){l+) [MeV] 

(£„r)(2+) [MeV] 

a,{\-) [fm] 

Q.(2-) [fm] 



TOSM 



Inert Core 



(0.22, 0.09) (0.03, 0.005) 

(0.64, 0.45) (0.33, 0.20) 

-5.6 1.4 

-17.4 0.8 



Table 8 The resonance energies E,- 
and the decay widths F of the d- 
wave resonance states in '"Li in unit 
of MeV. 



(£r,r)(i-) 
(Er,r){2-) 
{Er,r){3-) 
{Er,r)(A-) 



TOSM 



(5.84, 5.16) 

(5.81, 5.20) 

(6.57, 6.31) 

(5.30, 3.84) 



Table 9 5 used in the 'Li-« potential, and the energy differences AE between the (Ls)" and (Op)^ 
configurations of " Li in MeV. 



Inert core Pairing Tensor Present 



8 

AE 



0.066 
2.1 



0.143 
1.4 



0.1502 0.1745 
0.5 -0.1 



For the i-wave states, their scattering lengths Oj of the ^Li+n system show neg- 
ative values. In particular, the 2^ state shows a large negative value of Oj, which is 
comparable to that of the nn system (—18.5 fm) with the '^o channel |8J, and indi- 
cates the existence of the virtual 5-state near the ^Vl+n threshold energy. Therefore 
the inversion phenomenon in '"Li is reasonably explained in the present model. The 
order of 2^ and 1 ^ also comes from the attractive ^E component in the pn interac- 
tion. 

For comparison, we calculate '"Li without the core excitations of ^Li ("Inert 
Core"), namely, we adopt only the single Qp-Qli configuration for ^Li without the 
Pauli blocking effect explained in Fig.|2T|(d). In this case, we adjust the 5 parame- 
ter of the potential strength (1 + 5)V„, to be 0.066. In Table|2l the p-wave resonances 
are obtained just above the ^Li-nn threshold energy, and a^ values show small posi- 
tive values for both 1^ and 2^ states, which means that the virtual s- states are not 
located near the '^Li-nn threshold, and the s-p shell gap is large. These results mean 
that the Pauli-blocking nicely describes the spectroscopic properties of '"Li. 

The d-wave resonance states of '"Li are also predicted using the TOSM case of 
^Li, as shown in Table |8] whose excitation energies are higher than those of the p- 
states. The whole spectra of '"Li is summarized in Fig.[22]in comparison with the 
experimental data II103I . 
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Fig. 22 '"Li spectrum using 
TOSM and Inert core for 
*Li description 1591 . Exper- 
imental data is taken from 
Ref {m\. 
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5.6 ^^Li 

We solve "Li in a coupled ^Li+«+n model and show the detailed properties of "Li 
in Table [TOl the partial wave components P{{nlj)^) for halo neutrons, the various 
r.m.s. radius, the relative distance between halo neutrons {R„-„) and the core-2n dis- 
tance {Rc-2n), and the expectation value of the opening angle between two neutrons 
9 measured from the ^Li core, respectively. The case using TOSM for ^Li gives a 
large P{{ls)^) value, comparable to ^((pi/2)^) for halo neutrons and a large mat- 
ter radius R„j for "Li, which are enough to explain the observations. The case of 
"Inert Core" gives small ^((1.?)^) and small R,,, values, which disagree with the ex- 
periments. From the difference between two models, it is found that the tensor and 



60 


1 1 1 1 


^ 50 


■ 


[ '% - 


55 


( 


^ (^' Rm ^ 


■ 




■:j 40 


~ 


. P(s') P^ 




(b)l 


. ~ 


"c 30 


/ 




L 


J - 


\P 20 


■ / 








""" 


CO 


^^•~~-....,,^^^ / 








■"^m 


°- 10 


^Mj 




P(S'^) 









n 






_ 





3.7 

3.5 
3.3 
3.1 
2.9 
2.7 
2.5 



Inert core Pairing Tensor Present 



DC 



Fig. 23 (Color online) (l.s)^ probability P{s^) and matter radius R„, of "Li with four models in 
comparison with the experiments ((a) [3], (b) [T] and (c) 1581 ). The scale of P{.s- ) (R„,) is right 
(left) hand side. 
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Table 10 Ground state properties of "Li with S2,, = 0.31 MeV, wtiere two kinds of the '^Li de- 
scriptions of TOSM and Inert Core are shown, respectively. Details are described in the text. 





TOSM 


Inert Core 


Expt. 


^((pl/2)')) [%] 


42.7 


90.6 


— 


(1^1/2)' 


46.9 


4.3 


45±10(3] 


{P3/2)' 


2.5 


0.8 


— 


{dill? 


1.9 


1.3 


— 


{dsnf 


4.1 


2.1 


— 


if5/2? 


0.5 


0.2 


— 


[finf 


0.6 


0.3 


— 


Rm [fm] 


3.41 


2.99 


3.12±0.16(I1 
3.53±0.06|58| 
3.71±0.20|I04| 


Rp 


2.34 


2.24 


2.88±0.11(T] 


Rn 


3.73 


3.23 


3.21±0.17[a 


Rch 


2.44 


2.34 


2.467±0.037|105| 
2.423±0.034|106| 


Rn-n 


7.33 


6.43 




Rc-2n 


5.69 


4.26 




6 [deg.l 


65.3 


73.1 





pairing correlations in ^Li play important roles to break the magicity and make the 
halo structure of "Li, in addition to the s-p inversion phenomenon in '"Li. As was 
shown in the previous study [75], the Pauli blocking effect from the tensor correla- 
tion is stronger than the pairing one to break the magicity of "Li. 

In Fig. |23] "Present" is found to give a large amount of the (li)^ probability 
P{s^), 46.9% for the last two neutrons and a large matter radius R,n, 3.41 fm for "Li, 
which are enough to explain the observations. The probabilities of {pi/2)^^ {Pi/iY^ 
(c/5/2)^ and (t/3/2)^ for the last two neutrons are obtained as 42.7%, 2.5%, 4.1% and 
1 .9%, respectively. In Fig.|23] when we individually consider the tensor and pairing 
correlations for ^Li, P{s'^) is larger for the tensor case than for the pairing case. 
Finally, both blocking effects furthermore enhance P{s'^) and provide almost equal 
amount of (li)^ and {Qp)^ configurations. Hence, two correlations play important 
roles to break the magicity and make the halo structure for "Li. 

In Table|9] we also estimate the relative energy difference AE between (L?)" and 
(Op)^ configurations for "Li using the mixing probabilities of these configurations 
and the coupling matrix element between them as 0.5 MeV obtained in Ref. ||24| . 
The present model is found to give the degenerated energies enough to cause a 
large coupling between the (Op)^ and {Is)^ configurations by the pairing interaction 
between the last neutrons. 

In addition to the matter radius, the halo structure also affects the proton radius of 
"Li, because of the recoil effect of the cm. motion. In the three-body model of "Li, 
its proton radius {Rp) consisting of the proton radius of ^Li and the relative distance 
between ''Li and the cm. of two neutrons (/?c-2n) with the following relation 
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(^C^Li)) = {rICLi)) + i^-j (/??_2„), (86) 

where the second term represents the recoil effect. When the halo structure devel- 
ops, {R'^_2n) is expected to be large. Experimentally, considering the nucleon radius, 
the charge radius of "Li was measured recently and its value is 2.467±0.037 fm, 
which is enhanced from the one of '^Li, 2.217±0.035 fm II105I . The improved cal- 
culation for the isotope shift determination II 1 061 shows that 2.423±0.037 fm and 
2.185±0.033 fm for "Li and ^Li, respectively. The present wave functions provide 
2.44 fm and 2.23 fm for "Li and ^Li, respectively, which are in a good agreement 
with the experimental values. This enhancement is mainly caused by the large value 



of \/ {Rl_2n) obtained as 5.69 fm. For comparison, the distance between last two 
neutrons is 7.33 fm, which is larger than the core-2n case. 



5.7 Electromagnetic properties of Li isotopes 

We show the Q and jX moments of ^Li and "Li in Tables [TTI and [T2l where only the 
absolute value of the Q moment of "Li is reported in the experiments II107I . The 
present model describes reasonably those values for ^Li and "Li. For Q moments, 
the values of ^Li and " Li do not differ so much to each other This result is similar to 
that of the anti-symmetrized molecular dynamics (AMD) model 110811 and different 
from the stochastic variational method (SVM) based on the multi-cluster model 
lfT2ll . Here, similar to the charge radius, we discuss the recoil effect in the Q moment 
of " Li by expanding its operator Q ( "Li) into the core part (Q (^Li)), the last neutron 
part and their coupling part as 



^llT4^ _ ^^9t - /16?r, / 2 



e("Li) = QCU) + ^^3e i^-^j ^2o(Rc-2„) 

- Sn^^lOiCu),&^iiRc-2n)]20, (87) 

Oi,nCLi)^e ^ ^i^(a,), (88) 

/eproton 

where ^,„(r) = r'Yi„{r) and {a,} are the internal coordinates of protons in ^Li. 
In our wave function of "Li, last two neutrons almost form the 0^ state with the 
probability of around 99%. In that case, the Q moment for the relative motion part 
of ^Li-2«(c-2«) having the relative coordinate R^ -2n almost vanishes because of the 
non-zero rank properties of the Q moment operator. This means that the recoil effect 
from the clusterization is negligible. Therefore, the Q moment of "Li is caused 
mainly by the ^Li core part inside "Li. In our wave function, the spatial properties 
of the proton part of '^Li inside "Li do not change so much. Hence, the Q moment 
of "Li is similar to the value of the isolated ^Li. Small enhancement from ^Li to 
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' ' Li mainly comes from the lacking of the high momentum component of the tensor 
correlation due to the Pauli-blocking in ' 'Li, which extends the radial wave function 
of "Li. The experimental information of the Q moment is important to understand 
the structure of "Li. It is highly desired that further experimental data are available 
for the Q moment of "Li. 

For the jU moment, the value observed in "Li is almost the Schmidt value of 3.79 
jXn of the 0^3/2 proton. In ^Li, the p\i2 proton is slightly mixed, which decreases 
the /i moment. In ^Li, this pi/2 proton is excited from the O/93/2-orbit in a pair 
with the pij2 neutron and the 2p-2h excitation is Pauli -blocked in "Li due to the 
additional neutrons. As a result, the excitation oi pin is suppressed, which makes 
the jX moment of "Li close to the Schmidt value of the J33/2 orbit. The tendency 
of the increase of the ji moment from ^Li to "Li can be obtained also in the shell 
model analysis using various effective interactions IIIIOI . 

We further calculate the three-body Coulomb breakup strength of "Li into the 
^Li+n+n system to investigate the properties of the dipole excited states and com- 
pare the strength with the new data from the RIKEN group |16J. We use the Green's 
function method combined with the complex scaling method B3l to calculate the 
three-body breakup strength j28:| using the dipole strength and the equivalent photon 
method, where the experimental energy resolution is taken into account |6]. We do 
not find any resonances with a sharp decay width enough to make a resonance struc- 
ture. In Fig. |24] it is found that the present model well reproduces the experiment, 
in particular, for low energy enhancement and its magnitude. Seeing more closely, 
however, our results seem to underestimate the cross section at £ > 1 MeV, while 
overestimate at low energy peak region slightly. As a result, the integrated dipole 
strength forE < 3 MeV gives L35 e^fm", which agrees with the experimental value 
of L42 ±0.18 e^fm^lIS). We also slightly change the two neutron separation energy 
of the ground state of "Li, which is close to the recent observation 11111111121 . 

Summarizing this section, we have considered newly the tensor correlation in 
"Li based on the extended three-body model. We have found that the tensor and 
pairing correlations play important roles in ^Li with different spatial characteris- 
tics, where the tensor correlation prefers a shrunk spatial extension. The tensor and 
pairing correlations in ^Li inside "Li are then Pauli -blocked by additional two neu- 
trons, which makes the {Is)^ and (Op)^ configurations close to each other and hence 
activates the pairing interaction to mix about equal amount of two configurations. 



Table 11 Q moments of 'Li and " Li in units of e 
fm^. 



Table 12 fi moments of 'Li and '^Li 
in units of ^^r. 





"Li 


"Li 


TOSM 


-2.65 


-2.80 


AMDI108I 


-2.66 


-2.94 


svMfni 


-2.74 


-3.71 


ExDt.11071 


-2.74 ±0.10 


3.12±0.45(|e|) 


Expt.11091 


-3.06 ±0.02 


— 





'Li 


"Li 


TOSM 


3.69 


3.77 


AMDI108I 


3.42 


3.76 


SVMfni 


3.43 


3.23 


Expt. 11071 


3.44 


3.67 



Di-neutron clustering and deuteron-like tensor correlation in "Li 



55 



3.5 




Present (S2n=0.31MeV) — 
(S2n=0.38MeV) 



I®A ;.<" 



h<i> ,o^$^q'.<D^f$'^<. 



12 3 4 

Energy [MeV] 

Fig. 24 Calculated Coulomb breakup cross section of "Li into '^Li+n+n measured from the 
^Li+n+n threshold energy. Data are taken from Ref. (6). 



As a result we naturally explain the breaking of magicity and the halo formation 
for "Li. We also reproduce the recent results of the Coulomb breakup strength and 
the charge radius of '*Li. For '"Li, the inversion phenomenon is explained from the 
Pauli-blocking effect. 



6 Conclusion 



We have presented the physics of the di-neutron clustering and the deuteron-like 
tensor correlations by focusing on the halo structure of "Li. The halo structure pro- 
vides an ideal platform for the di-neutron clustering to play an important role. This 
di-neutron clustering phenomenon is strongly related with the central interaction in 
the ' So channel, where the large scattering length in the nucleon-nucleon scattering 
suggests the appearance of the bound state by changing slightly the environment. 
On the other hand, in order to produce the halo structure, there should be an active 
participation of the Sin configuration in the neutron wave function. This participa- 
tion of the 51/2 orbit was very difficult in the standard shell model framework. We 
had to invoke the deuteron-like tensor correlation in "Li, which blocks the two neu- 
trons to enter the pi^ orbit and hence provides a mechanism to put neutrons in the 



h/2 



orbit. 



The theoretical description of these two new correlations to realize in "Li was 
very difficult. In fact, we had to develop the cluster orbit shell model (COSM) to 
handle an extended object and further the hybrid-Vr model (Hybrid-V T) to treat the 
di-neutron clustering correlation in finite nuclei. We had to deal with unbound and 
resonance states quantitatively, since the halo structure appears when the binding 
energy of the last neutrons is very small. The excitation function of several MeV 
forces us to treat the spectral function in the continuum region. At the same time, the 
neighboring nuclei are also unbound. Hence, we had to develop the complex scaling 
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method (CSM) to treat the continuum and resonance states. All these theoretically 
involved methods have been developed for the quantitative description of the halo 
nucleus "Li. The most essential mechanism to bring down the ij « configuration 
for the halo formation came from the deuteron-like tensor correlation in nuclei. The 
strong binding energy of ^He comes from the strong tensor interaction. We had 
to develop the tensor optimized shell model (TOSM) to describe the deuteron-like 
tensor correlation in the shell model basis. The TOSM provided a clear key to push 
up the configuration involving the p^n orbit, since the strong binding of the ^He 
nucleus requires the use of the pi/2 orbit. Hence, naturally the ,vi/2 configuration is 
energetically favored for the halo structure formation. 

We are used to treat the central interaction with the spin-orbit interaction in the 
shell model basis to treat many body systems. In the theoretical challenge to de- 
scribe the halo structure of "Li, we had to face to treat the strong tensor interaction 
for the important role of the deuteron-like tensor correlation in many body systems. 
We have seen that the deuteron-like correlation worked out in the shell model basis 
is able to treat the Pauli blocking effect due to the {pinY neutrons to wash away 
theN=8 shell gap. This effect allows the participation of the (^1/2)^ configuration to 
provide a platform to develop the di-neutron clustering correlation in the halo struc- 
ture of "Li. Hence, we expect many interesting many-body phenomena in unstable 
nuclei to be found in near future, where the deuteron-like tensor and/or di-neutron 
clustering correlations play important roles. 

In this lecture note, we have gone through all these theoretical materials in de- 
tails. We have tried to make the motivation of the development of the theoretical 
tools in each step to describe the halo nucleus quantitatively. These theoretical 
frameworks are not only for use of the halo structure, but should play a very im- 
portant role for the description of finite nuclei. In fact, the TOSM is essential to 
describe the deuteron-like tensor correlation, which should be the most important 
ingredient to provide large binding energies for all the nuclei. 
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